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ABSTRACT 

Context. EX Lupi is a young, accreting MO star, prototype of EXor variable stars. Its spectrum is very rich in emission lines, including 
many metallic lines with narrow and broad components. It has been also proposed to have a close companion. 

Aims. We use the metallic emission lines to study the accretion structures and to test the companion hypothesis. 

Methods. We analyse 54 spectra obtained during 5 years of quiescence time. We study the line profile variability and the radial 
velocity of the metallic emission lines. We use the velocity signatures of different species with various excitation conditions and their 
time dependency to track the dynamics associated to accretion. 

Results. We observe periodic velocity variations in the broad and the narrow line components, consistent with rotational modulation. 
The modulation is stronger for lines with higher excitation potentials, which are likely produced in a confined area, very close to the 
accretion shock. 

Conclusions. We propose that the narrow line components are produced in the post-shock region, while the broad components 
originate in the more extended, pre-shock material. All the emission lines suffer velocity modulation due to the rotation of the star. 
The broad components are responsible for the line-dependent veiling observed in EX Lupi. We demonstrate that a rotationally- 
modulated line-dependent veiling can explain the radial velocity signature of the photospheric absorption lines, making the close-in 
companion hypothesis unnecessary. The accretion structure is locked to the star and very stable during the 5 years of observations. Not 
all stars with similar spectral types and accretion rates show the same metallic emission lines, which could be related to differences in 
temperature and density in their accretion structure(s). The contamination of photospheric signatures by accretion-related processes 
can be turned into a very useful tool to determine the innermost details of the accretion channels in the proximities of the star. 
The presence of emission lines from very stable accretion columns will nevertheless be a very strong limitation for the detection 
of companions by radial velocity in young stars, given the similarity of the accretion-related signatures with those produced by a 
companion. 

Key words, stars: pre-main sequence, stars: variability, stars: EX Lupi, protoplanetary disks, accretion, techniques: spectroscopic, 
radial velocities 


1. Introduction 


EX Lupi is the prototype of EXor variables and a remarkable ob¬ 
ject. Its strong variability episodes are known for more than 50 
years (McLaughlin 1946; Herbig 1950). The photospheric emis¬ 
sion observed during the quiescent phases is compatible with a 
young star with spectral type MO (Herbig 1977; Gras-Velazquez 
& Ray 2005). The star is known to have short-timescale, mod¬ 
erate (1-2 mags) variability (Herbig 1977; Lehman et al. 1995) 
caused by small variations in the accretion rate, together with 
more rare extreme variability episodes (increasing by 4-5 mags 
in the optical), of which only two have been documented to date: 
the first one in the 50’s (Herbig 1977), and a second one in 2008 
(Jones 2008). In these extreme variability episodes, the accre¬ 
tion rate increases by more than two orders of magnitude dur¬ 


ing a few months time, before the star goes back to quiescence 
(Abraham et al. 2009). 

The main characteristic of the EX Lupi spectra is the large 
number of emission lines observed, both in quiescence and in 
outburst (Patten et al. 1994; Herbig et al. 2001; Sipos et al. 2009; 
Kospal et al. 2008, 2011; Sicilia-Aguilar et al. 2012, from now 
on SA12). During its strong 2008 outburst, the whole spectrum 
was dominated by emission lines and the stellar photosphere was 
completely veiled, with the Li I 6708A line being the only one 
that still showed a recognisable absorption (SA12). Such rich 
emission-line spectra are very rare among classical T Tauri stars 
(CTTS) and even other EXors. To our knowledge, the only other 
object that shows a similar outburst spectrum is the M5-type 
EXor ASASN13db (Holoien et al. 2014). The number of identi¬ 
fied emission lines in the optical is of the order of a thousand, al¬ 
though there are many more lines that are hard to identify due to 
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blends (SA12). Most of the emission lines show broad and nar¬ 
row components (BC and NC), with the remarkable exception of 
the hydrogen lines, which have only broad, very complex, pro¬ 
files. The BC shows a strong day-to-day modulation, consistent 
with bulk motions of clumpy, non-axisymmetric material rotat¬ 
ing and accreting onto the star (SA12). In quiescence, the num¬ 
ber of emission lines is reduced to about 200, although some 
photospheric lines show cores consistent with weak, NC emis¬ 
sion (SA12). Most of the metallic lines in quiescence are narrow, 
with very weak or absent BC, although the typical CTTS emis¬ 
sion lines are broad (e.g. the hydrogen Balmer series) or show 
narrow and broad components (e.g. Ca II, some of the strongest 
Fe II lines, the He I 5875A line). Except for the Ca IIIR triplet, 
most of the BC of metallic and He I lines are weak during qui¬ 
escence. All the observed emission lines correspond to permit¬ 
ted transitions, with no evidence for forbidden emission lines (in 
both outburst and quiescence). 

The photospheric absorption lines during quiescence show 
a strong radial velocity (RV) modulation, which has been inter¬ 
preted as caused by a brown dwarf (BD) companion, although 
hot-spot-induced RV signatures could not be excluded (Kospal 
et al. 2014; from now on K14). There are several precedents 
of objects where the RV suggested initially the presence of a 
companion, which was later on ruled out as a signal caused by 
accretion and/or strong activity, such as RW Auriga (Gahm et 
al. 1999; Petrov et al. 2001; Dodin et al. 2012) and RU Lupi 
(Stempels et al. 2007; Gahm et al. 2013). In fact, using radial 
velocity techniques to detect companions to young (active, vari¬ 
able, accreting) stars has well-known problems (Crockett et al. 
2012; Jeffers et al. 2014; Dumusque et al. 2014). 

The emission lines observed in EX Lupi are a key to dis¬ 
entangle accretion and companion effects and to understand the 
way accretion proceeds in the EXor prototype. Emission lines 
have been used for more than 30 years to study the circumstellar 
environment, accretion and winds in CTTS, specially consider¬ 
ing the strong Ha, He I, and Na I D lines, with the work started 
by Hartmann (1982), Appenzeller et al. (1983), Edwards et al. 
(1994), Hamann & Persson (1992), and Hartmann et al. (1994) 
and continued by many others. In contrast with the more com¬ 
plex strong lines, the weak metallic lines (Ee I/II, Ti II), have less 
complex profiles, not affected by self-absorption. The metallic 
lines also cover a broad range of excitation conditions, with ex¬ 
citation potentials ranging from few eV to several tens of eV and 
different transition probabilities. This gives us a unique chance to 
probe regions of the accreting system that show different values 
of temperature and density, constructing a 3-dimensional pic¬ 
ture of the EX Lupi environment. The day-to-day modulation 
observed in the broad emission lines during the 2008 outburst 
showed clear velocity differences depending on the species, e.g. 
comparing ionised vs neutral lines, or lines with high excitation 
potentials such as He I vs lines with weaker excitation poten¬ 
tials such as Ca II or Ee I (SA12). The differences between lines 
were consistent with the differences in density and temperature 
expected along an extended (~0.1 AU) accretion column(s) or 
non-axisymmetric rotating/infalling structure(s) (SA12). A sim¬ 
ilar approach, but using absorption lines, has been used to track 
cold matter moving around intermediate-mass stars (Mora et al. 
2002, 2004), and to trace the physical conditions along the ac¬ 
cretion column in S CrA SE (Petrov et al. 2014). With the time 
coverage we now have for EX Lupi in quiescence, we can not 
only explore the physical properties, but also use the time vari¬ 
able to explore the dynamics of the accretion columns. 

In this paper, we analyse the emission lines of the quiescent 
EX Lupi, turning the classical problem of accretion-related sig¬ 


natures affecting the radial velocity measurements into a tool to 
investigate the structure of its accretion columns. Section [^in¬ 
troduces the EX Lupi spectroscopy data. In Section]^ we ex¬ 
plore and quantify the dynamics observed in narrow and broad 
emission lines. Simple dynamical and physical models for the 
accreting structure are presented in Section]^ Section [^brings 
together the observations within a unified picture of rotational 
modulation. Einally, Sectionj^summarises our results. 

2. Observations, data reduction, and disk and 
steiiar properties 

2.1. Observations and data reduction 

EX Lupi was monitored using the ’’Eiber-fed Extended Range 
Optical Spectrograph” (EEROS; Kaufer et al. 1999) mounted 
on the 2.2m telescope and the ’’High Accuracy Radial veloc¬ 
ity Planet Searcher” (HARPS; Mayor et al. 2003) spectrograph 
on the 3.6m telescope, both in La Silla, Chile. The observations 
were part of various observing programs (see K14 for details on 
the quiescence data). The spectroscopic followup of EX Lupi 
consist of 54 spectra (see Table [AT] in Appendix [Aji taken dur¬ 
ing the quiescence phase of EX Lupi. These spectra correspond 
to the high S/N observations listed by K14. The observations 
cover about 5 years of quiescence data taken between 2007 and 
2012 (excluding the data taken during the 2008 EX Lupi out¬ 
burst; discussed in SA12). The data were taken at irregular inter¬ 
vals, ranging from two spectra per night to few days/months in¬ 
tervals, depending on telescope availability. The exposure times 
ranged between 10-50 minutes, depending on the brightness of 
the star and the atmospheric conditions. EEROS has a resolu¬ 
tion of R=48000, and a wavelength coverage between 3500- 
9200A, containing in total 39 echelle orders with some small 
gaps. HARPS data have a higher resolution (R=l 15000) but a 
more reduced spectral coverage (3780-6910A). Eor the line com¬ 
parison in this study, the HARPS spectra have been resampled to 
the EEROS resolution, which also allows us to reach a more uni¬ 
form signal-to-noise ratio (S/N). 

The data were reduced using the automated HARPS and 
EEROS pipelines, and then corrected to remove the offset in the 
barycentric velocity for the EEROS specra (Muller et al. 2013) 
and to remove the radial velocity instrumental offset betwen the 
EEROS and the HARPS spectra of -0.291 km/s (K14). Both ef¬ 
fects are very small and thus have little effect on the broad emis¬ 
sion lines, but they could have some influence on the radial ve¬ 
locity measurements of the narrow emission line components, of 
the order of few km/s. 

2.2. Steiiar and disk properties revisited 

The spectra of EX Lupi in quiescence are consistent with a MO 
star (stellar mass ~0.6 Mq; Gras-Velazquez & Ray 2005) with 
numerous emission lines. The luminosity is highly variable due 
to variations in the accretion rate. The estimated stellar radius 
is ~1.6 Rq (Sipos et al. 2009). The rotational velocity (ysin/) is 
consistent with the value of 4.4+2.0 km/s (Sipos et al. 2009). 

The lack of near-IR excess is interpreted as a small, dust- 
depleted hole in the disk, with size ~0.3AU (Sipos et al. 2009). 
CO observations (Goto et al. 2011; Banzatti et al. 2015) are also 
in agreement with a sub-AU hole, mostly devoid of dust. The 
disk had been suggested to be relatively massive (0.025 Mq for 
a star with mass ~0.6Mq) based on disk models containing only 
small dust grains (Sipos et al. 2009). If instead we assume a more 
typical maximum grain size of 100 pm as it is usually found 
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in other T Tauri stars (D’Alessio et al. 2001, 2006; Andrews & 
Williams 2007), the disk mass does not need to be so extreme 
and could be below 1% of the stellar mass (Appendix [B] i. 

The typical accretion rate during the quiescence periods has 
been estimated to be l-3xl0“*°Mo/yr (SA12, based on Ha 10% 
widths; Natta et al. 2004). Estimates of the accretion rate based 
on Pa j3 (Sipos et al. 2009) and on multiple emission lines (us¬ 
ing the relations of Alcala et al. 2014; Appendix [C) result in a 
value of 4+2xlO“^°Mo/yr. The observed Ha width at 10% of the 
peak (Natta et al. 2004) during our quiescence data are in agree¬ 
ment with these values, although for some of the spectra with 
stronger lines, the accretion rate increases up to ~10“®Mo/yr. 
The data discussed here thus includes variations in the accretion 
rate within one order of magnitude, well below the outburst case 
where the accretion rate increases by more than two orders of 
magnitude (SA12). 

3. Quantifying the iine variabiiity for a dynamicai 
anaiysis 

The emission lines are very complex. To hnd a proxy to measure 
line variability and dynamics, we take two different approaches: 

- First, we explore the pixel-by-pixel variation of the broad 
components (BC) without assuming any particular model for 
the line. This method was used by Lago & Gameiro (1998) 
and Alencar et al. (2001) to study variability patterns in T 
Tauri stars. It allows us to identify the dynamical effects that 
dominate the changes in the line profiles, which parts of the 
line have the highest variability, and whether variability af¬ 
fects all parts or velocity components of the line similarly. It 
also allows us to check for correlations between the different 
velocities and lines. 

- Second, we use a multi-Gaussian fit to reproduce the line 
emission and explore the variations in intensity, radial veloc¬ 
ity, and width, of the various Gaussian components. This is 
particularly useful to study the behaviour of the narrow com¬ 
ponents (NC), whose velocity, amplitude, and width do not 
depend on the particular choice of Gaussian ht. Since most 
of the metallic lines have weak (often undetectable) BC with 
poorly resolved prohles, the metallic line analysis refers al¬ 
most exclusively to the NC of the lines. 

3.1. A pixel-by-pixel dynamical view of the broad component 

To explore the pixel-by-pixel variation of the emission lines, we 
hrst resampled all the spectra to the (lower) EEROS resolution. 
The richness of features in the EX Lupi spectra is such, that 
hnding a proper normalisation at all wavelengths is very hard. 
Therefore, for each line, the spectra are normalised to the local 
continuum levels, measured on both sides of the line on parts 
of the spectrum that are not affected by emission or absorption 
lines. We then estimated the median and standard deviation of 
the emission line at each pixel, after clipping to remove bad pix¬ 
els and cosmic rays. The procedure was repeated for all the broad 
lines along the extent of the BC line wings. 

We also explore to which extent the emission at a given ve¬ 
locity is correlated with the emission at any other velocity within 
the same line, and how different lines are correlated with each 
other. If we attribute each velocity bin to a particular gas par¬ 
cel (or parcels) moving at the same rate, the autocorrelation ma¬ 
trix and the two-line cross-correlations reveal which parts of the 
system may be connected. Nevertheless, there could be material 
moving at the same projected velocity, but coming from different 


physical locations, which can alter the pixel-by-pixel velocity 
correlations. In addition, Dupree et al. (2012) demonstrated that 
there are signihcant temporal variations of the order of hours be¬ 
tween the various accretion indicators (including the B aimer H 
emission, the veiling continuum, and the He I lines) in a dumpy- 
accretion episode. Since the reaction time differences between 
various lines are longer than our exposure times, this could wash 
out the correlations. 

Figures and show both the line variance and median 
value and the autocorrelation matrix at different velocities for 
hydrogen and Ca II lines. The Hydrogen Balmer lines are, as ex¬ 
pected (Edwards et al. 1987; Appenzeller et al. 1989; Muzerolle 
et al. 2001), very complex and asymmetric. They often have 
self-absorptions that may correspond to colder wind/infall com¬ 
ponents, which can dilute the correlations or otherwise produce 
very complex patterns. Ha is particularly remarkable, with a pro¬ 
file with both blueshifted and redshifted absorptions. The blue 
part of the spectrum shows a higher variability, probably due to 
the presence of a variable wind. Part of the redshifted emission 
may also be suppressed by self-absorption in an optically thick 
accretion flow. The maximum velocities observed in the red and 
in the blue increase when the strength of the line increases, as ob¬ 
served during outburst and expected in the case of fluctuations 
of the accretion rate (Natta et al. 2004). If attributed to free-fall 
(infall-dominated) alone, the maximum velocities (+300 km/s) 
would suggest material flowing onto the star from distances of 
~2 stellar radii, considering Mt=0.6Mo and Rt = 1.6Ro. There is 
also a tendency of the line to become more symmetric when its 
intensity (and thus accretion rate) decreases, also a sign of the 
complex profile being strongly affected by self-absorption. The 
higher Balmer lines become increasingly symmetric, with sim¬ 
ilar variability on their blue- and redshifted sides. For the Ha 
and Hy6 lines, we hnd a strong correlation between the blue- and 
red-shifted parts, while the correlation of the line wings with the 
central part of the line is signihcantly lower. This is a further sign 
that the zero-velocity components may be saturated and/or pro¬ 
duced in several spatially different or very extended locations. 

The Ca II lines, especially the IR triplet (only covered by 
EEROS) also show very broad components. The Ca II IR triplet 
BCs have a remarkable prohle reminiscent of the outburst pro¬ 
hles of the metallic lines (SA12): the BC of the line changes 
from blue- to redshifted within timescales of days, spanning a 
large range of velocities. The BC also decreases to be nearly un¬ 
detectable when the accretion rate decreases (see Appendix 
The median line prohle is, in strong contrast to the individual 
lines, very symmetric, as well as the pixel-by-pixel standard de¬ 
viation, showing that the variability in the blue and in the red is 
very similar, but independent. The Ca II 8542A line falls par¬ 
tially in the EEROS spectrograph gap, but it appears similar the 
other two Ca II IR triplet lines. The Ca II IR triplet is thus con¬ 
sistent with bulk motions of parcels of gas that are sometimes 
blueshifted and sometimes redshifted, rather than to related pro¬ 
cesses (such as infall/accretion and accretion-powered wind) that 
produce correlated blueshifted and redshifted features. At hrst 
order, the lack of anticorrelations is also a problem for the com¬ 
panion scenario, where we would expect the RV shift to increase 
one side of the line while decreasing the other one, although the 
shift could be signihcantly smaller than the BC velocity widths 
and thus hard to detect. Although the Ca II line forms closer to 
the star, this asymmetric motion is consistent with the CO obser¬ 
vations of Goto et al. (2011) and Banzatti et al. (2015) and with 
the day-to-day variations in the line prohles during the 2008 out¬ 
burst (SA12), suggesting that the asymmetric, rotating/infalling 
structure that hlls in the inner disk probably continues onto the 
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Fig. 1. Ha and H/l line profile variability (left) and autocorrelation matrix (right). The pixel-by-pixel median of the line is shown 
as a blue thick line. The pixel-by-pixel rms is shown as the thick black line and shaded area. The zero velocities/wavelengths of 
the lines are marked as dashed vertical lines. The correlation coefficient (estimated as a Spearman rank correlation with value r and 
false alarm probability p) is shown on the right in the color scheme, with the black contours marking the high significance (p>10“^) 
parts with 1-0.6,0.7,0.8,0.9. 


inner magnetospheric accretion region. The Ca II H and K lines 
have asymmetric profiles with a blueshifted absorption that is in 
part correlated with the wings on the redshifted side of the line, 
although the velocity range is small compared to the H B aimer 
lines or the Ca IIIR tripIejM 

We can take a further step exploring the connections between 
different lines and velocities by calculating the pixel-by-pixel 
cross-correlation between different lines (Figure [^. The cross¬ 
correlation matrix for the Ca II 8498 and 8662A shows the same 
pattern as the autocorrelation matrix of each line, a sign that 
they are nearly identical, as expected. There is an anticorrela¬ 
tion between Ca II K and Ca II IR, likely due to S/N effects and 
variations of the accretion rate. A higher accretion rate produces 
strong Ca II IR, but also a higher blue continuum, and thus a 
lower peak over the continuum of all the UV and blue lines. At 
low accretion rates, the excess blue continuum decreases, and 
since EX Lupi is a MO star with little blue/UV emission, the con¬ 
tinuum at short wavelengths is noise-dominated and the strength 
of blue/UV lines over the continuum is artificially increased. A 
similar anticorrelation appears between Ha and Call K: the Ha 
blue and red wings are anticorrelated with the Ca IIK line. Since 
the red and blue Ha wings (and the near-UV continuum) grow 


' Note that the Ca II H line is contaminated by He. 


when the accretion rate increases, the Ca II K peak over the con¬ 
tinuum decreases if the Ca II K line flux does not increase at 
the same rate. This is a sign that we need to be careful when 
comparing lines at different wavelengths. 

Cross-correlation reveals that the wings of Ha and the BC 
of the Ca II IR triplet lines are correlated (Figure]^, a sign that 
the strength of the Ca II IR triplet BC depends on the accretion 
rate, although the velocities and line profiles change in an inde¬ 
pendent way. There is no significant correlation between the Ha 
peak and the Ca II IR lines, probably because Ha is strongly sat¬ 
urated. The cross-correlation between the Ca II IR triplet and the 
He I lines does not give any significant result. This suggests that 
both lines respond to different physical processes that are likely 
to happen in different locations and at slightly different times 
(Dupree et al. 2012). 

3.2. Multi-Gaussian fit and radial velocities 

The second step involves studying and quantifying the strength, 
width, and radial velocity of the different emission lines. We fit¬ 
ted small sections of the normalised and continuum-subtracted 
spectrum (typically, +100 km/s around the line for narrow fea¬ 
tures, and up to +300 km/s around features with BC) with a 
multi-Gaussian model. For the strongest, complex lines (Ha, H/1, 
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Fig. 2. Same as Figure [^for the Ca II K (top) and 8498A (bottom) lines. The pixel-by-pixel median of the line is shown as a blue 
thick line. The pixel-by-pixel rms is shown as the thick black line and shaded area. The zero velocities/wavelengths of the lines are 
marked as dashed vertical lines. The correlation coefficient (estimated as a Spearman rank correlation with value r and false alarm 
probability p) is shown on the right in the color scheme, with the black contours marking the high signihcance (p>10“^) parts with 


r=0.6,0.7,0.8,0.9. 


etc), any attempt of simple multi-Gaussian htting or interpreta¬ 
tion is impossible. The He 1 and He II lines are less complex, 
as well as the metallic lines (Ca II, Fe I, Fe II, Ti II, Si II, Mg I, 
among others). These weaker emission lines are the main subject 
of this section. 

Our aim is not to hnd the best mathematical ht, but a simple 
and robust one that can be used for all the data taken at differ¬ 
ent dates in order to compare them. As a result, a direct physi¬ 
cal interpretation may be difficult, but we will be able to detect, 
quantify, and explore differences in velocities, widths, and am¬ 
plitudes in the spectral lines. A multi-Gaussian ht is strongly 
degenerated, so we imposed a few restrictions based on the 
visual aspect of the lines and continuum to make the models 
comparable. We used 3 Gaussians with independent amplitudes, 
widths, and wavelength offsets. In some cases where more than 
one line is present within the considered section, we added a 
4th Gaussian component (for instance, to ht the O I triplet at 
7774A). The lines have a NC-i-BC structure, so we request one 
of the Gaussians to be narrow, with a width between 8-30 km/s 
(or two narrow Gaussians in case of two nearby narrow lines 
such as Fe I 5269/70A). The Gaussian amplitudes are forced to 
be positive, except for the He I lines, which are better htted by in¬ 
cluding one Gaussian component in absorption, and cases where 
the emission line appears on top of a photospheric absorption 


feature (e.g. Fe I 5269/70A). All hts were visually inspected, re¬ 
moving any parts of the line or continuum that were affected by 
spikes. Figure|^shows some examples of the multi-Gaussian hts 
for various lines. 

We hnd that the NC ht is independent of the assumptions 
made on the htted Gaussians. Nevertheless, htting the BC is 
strongly degenerated, as most spectra can be equally well ht¬ 
ted by various very different Gaussian components. Therefore, 
the results of the BC htting have to be handled with care, while 
the amplitudes, velocities, and widths of the NC are very robust. 
Details on individual lines are given in the Appendix]^ 

3.2.1. Analysis of the narrow component 

Our detailed multi-Gaussian ht allows us to obtain more accurate 
velocity measurements than the automated hting and averaging 
method in K14, allowing us to treat each line independently. Our 
discussion here is thus focused on the velocity modulation of 
the strong NC emission lines. Weaker narrow emission lines are 
noisier, although they are consistent with the same modulation. 
K14 noted that the NC of the emission lines shows radial ve¬ 
locity patterns that are off-phase with respect to the RV curves 
of the photospheric absorption lines. To quantify the velocity 
modulation of the NC emission, we run a generalised Lomb- 
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Fig. 3. Pixel-by-pixel velocity cross-correlations between different lines. From top to bottom, left to right: Ca II 8498 vs 8662A, Ca 
II K vs Ca II 8498A, Ca E 8498A vs Ha, Ca E 8222A vs He I 5875A, Ha vs He I 5875A and Ca E 8498A vs H/3. The correlation 
coefficient (estimated as a Spearman rank correlation with value r and false alarm probability p) is shown on the right in the color 
scheme, with the contours marking the high significance (p>10“^) parts with r=0.6,0.7,0.8,0.9. 


Scargle periodogram (GLSP; Scargle 1982; Horne & Baliunas 
1986; Zechmeister & Kuerster 2009) to search for periodic sig¬ 
nals. Each line has a very stable characteristic line width, which 
can be as narrow as 9 km/s for some Ti E lines and as broad as 20 
km/s for He I lines. The data were clipped based on the NC line 
width, since significantly broader lines are usually contaminated 
by the BC, while very narrow ones are often affected by spikes 
or low S/N, leading in both cases to inaccurate NC velocities. 


The errors in the line fit can be very complex, considering that 
the sum of Gaussians is a degenerated fit. Therefore, we estimate 
the errors in the velocity as being proportional to the line width 
and inversely proportional to the line strength, including a factor 
to account for the finite spectral resolution, resulting in typical 
errors well below 0.5 km/s. The clipping procedure leaves us be¬ 
tween 30 to over 50 points per emission line, with the strongest 
lines having more valid points than the faint, low S/N ones. The 
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Fig. 4. Some of the multi-Gaussian fits for one of the dates (JD 2455674.82443). The normalised, continuum-subtracted observed 
spectra are displayed as thick black lines. Individual Gaussian components are plotted as thin lines (colours; green, red, blue, cyan). 
The combined fit is plotted as a magenta thick line. The zero velocities of the central line are also marked with a vertical dotted line. 
The plots are labelled with the line(s) elements and wavelength. Note that the very complex shape of the continuum around the Fe I 
5269/70A lines is better fitted by two strong Gaussians, one in emission and the other in absorption (partly out of the figure). These 
BC have no physical meaning/interpretation, but allow us to extract properly the strong NC of both lines. 



Period (d) 




Period (d) 


Period (d) 



Period (d) 


Fig. 5. Some examples of GLS periodograms for different lines. The 7.41d period (marked as a blue dashed line) becomes increas¬ 
ingly significant as we move to lines with higher excitation potentials. The FAP are marked as dotted horizontal lines and labelled 
accordingly. The situation observed for Mg I is typical of the rest of metallic lines labelled as having several peaks. 


results of the GLSP (periods and false-alarm probabilities, FAP) 
are summarised in Table IA.2I 

We find that some of the lines show significant periods sim¬ 
ilar to the radial velocity period of the photospheric absorption 
lines, 7.41d (K14). In most cases, rather than a single period, we 
observe a strong modulation and a signature of quasi-periodicity, 
with several similarly-significant peaks around the 7-8d range. 
Periods of 8.09 and 7.19d are frequent, and could be aliases of 
the 7.41d period. The strongest signature is observed for the He 
II 4686A line, for which we recover the period obtained for the 
photospheric absorption lines by K14, although this line is not 


particularly strong compared to others. Some of the strongest 
lines, such as the Ca II 8498 and 8662A lines, show no signifi¬ 
cant periodicity (as expected from the lack of modulation; Figure 
1^. Therefore, the presence (or lack of) significant periodic sig¬ 
natures is not an issue of S/hQ but rather due to the underlying 
structure of the emitting material. Figure shows some exam¬ 
ples of the resulting periodograms, and details about every indi¬ 
vidual line are summarized in the Appendix [A| Table 

^ Considering the consistent opacity-related shift observed between 
the two Ca IIIR lines, we should be able to detect any modulation with 
amplitude at the ~0.5 km/s level, if there was any. 


A.2 
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Fig. 6. Phase-folded curves (using the photospheric absorption 
line period and phase, which is nearly identical to the He II 
4686A one)) for several of the observed NC of lines. Note that 
two 7.417d periods are plotted for clarity. The dotted line rep¬ 
resents the RV curve measured for the photospheric absorption 
lines (K14). 


Plotting the NC central velocities against the He II (or K14) 
period reveals that most of the lines are consistent with similar 
sinusoidal periodic modulations with the same period than the 
RV signatures of the absorption lines (Figure |^. The RV mod¬ 
ulation depends on the line considered, ranging from lines that 
show hardly any modulation (such as the NC of the Ca II IR 
triplet), to others where the amplitude of the RV modulation is 
much larger than the RV amplitude of the photospheric absorp¬ 
tion lines (e.g. He I, He II show modulations up to 8 km/s peak- 
to-peak amplitude, compared to 4.36 km/s for the photospheric 
absorption lines). Among the parameters in the Gaussian fit (RV 
of the center of the Gaussian, amplitude, and width), only the 
RV shows a significant modulation. The line width is remark¬ 
ably stable within the same species and transition, and the am¬ 
plitude shows low-significance modulations that could be related 
to orientation of the emitting structure (Section]^. Interestingly, 
although the BC appear more intense during higher accretion 
phases, there is no evidence that the variations of the accretion 
rate for the quiescence data produce any significant change in 
the properties and velocities of the NC (see Appendix [C|). 

Most of the lines are redshifted by different amounts (from 
less than ~ 1 km/s for the Ca II IR lines, up to ~5 km/s for He 
lines) with respect to the zero-point stellar velocity (-0.52 km/s; 
K14). This suggests that the lines are not produced at the stel¬ 


lar photosphere, but in a slightly infalling structure. Since the 
redshifted velocities of the lines also vary from element to el¬ 
ement, the various species are not produced in the same place 
of the structure. Significantly redshifted He I lines (compared to 
the stellar velocity and to Fe I/II lines) have been also observed 
in other stars by Beristain et al. (1998, 2001) and Gahm et al. 
(2013). 

Figure]^ also reveals a systematic offset between lines from 
the same multiplet. This is particularly remarkable for the 8498 
and 8662A lines of the Ca II IR triplet, and can be explained by 
differences in the optical thickness of the lines. Although the 
third component of the triplet, at 8542A is very close to the 
FEROS spectrograph gap, the NC is often within the FEROS 
coverage, displaying a behaviour comparable to the 8498A line 
in peak intensity and zero-point velocity, while the 8662A line is 
redshifted by ~lkm/s with respect to them and has a lower peak 
intensity. The 8498A NC line has a tendency to being narrower 
than the other two. Hamann & Persson (1992) noticed the dif¬ 
ference in intensity and thickness of the 8498A vs 8542-8662A 
lines in CTTS (assuming a chromospheric origin for the lines) 
and offer two explanations: The 8662A reduced intensity could 
be caused by its upper level being partly depopulated by interac¬ 
tions between the Ca IIH line (which shares the same level) and 
the He line. Regarding the line thickness (and, in our case, also 
probably the offset velocity variation), an explanation may lie in 
the 9-times lower opacity of the 8498A line, compared to 8662A. 
This can cause a change in the line profile (the optically thicker 
line is broader) and also means that the 8498A line would be pro¬ 
duced in a deeper place, compared to the 8662A line. Hamann & 
Persson (1992) suggest a temperature inversion causing a source 
function inversion in the zone where the Ca II IR triplet orig¬ 
inates, which would mean that the 8498A line would be pro¬ 
duced in the deepest place, the 8662A line would be produced 
at the minimum, and the 8542A line would originate where the 
temperature raises again. Such a situation could also arise in a 
post-shock cooling region (Section |4.1| i. A similar explanation 
could also fit the differences between the Mg 15167/5172/5183A 
triplet, which also shows differences in strength and radial veloc¬ 
ity between the 5167A line and the rest, although the 5167A line 
coud be in part affected by a strong, nearby Ee I line. 

3.2.2. Analysis of the broad component 

We also searched for periodic signatures in the strong, rapidly 
variable BC of the Ca II 8498 and 8662A lines. The span of 
peak velocities is over ±100 km/s, and similarly good fits may 
show variations by several tens of km/s. We first compared the 
independent fits for the 8498 and 8662A lines to check how reli¬ 
able and stable the BC velocity determinations are. We estimated 
the Spearman rank correlation for all the three Gaussian fits 
components: amplitude, central wavelength, and width. Cross- 
correlating the results for the two lines (see Appendix [D| for de¬ 
tails) shows that there is a very strong correlation between the 
individual parameters of both lines. This means that we are ob¬ 
taining similar fits for these two very similar lines, despite vari¬ 
ations in the noise and continuum between both. Eor most of our 
data points, a single broad Gaussian dominates the BC fit, and 
the second broad Gaussian component is often marginal and/or 
poorly constrained, which leads to degeneracy. 

Eor a given line, we checked the cross-correlation between 
the fit parameters, to investigate how changes in one of the com¬ 
ponents (NC or the two broad Gaussian components) are re- 
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fleeted in the other ones. Appendix [D| details the procedure and 
results. We hnd that the NC is relatively independent in ampli¬ 
tude and central velocity from the BC. Only the width of the NC 
seems to increase slightly when the intensity of the BC increases. 
This could be a sign of an increase of the local turbulence at 
near-zero velocities for higher accretion rates (expected by com¬ 
parison with the outburst observations; SA12), but it could also 
reflect contamination of the NC by a more intense BC. The cen¬ 
tre and the amplitude of the two BC Gaussians are correlated, so 
their velocity increases when the amplitude of the BC increases. 
This is similar to what we observed during outburst: at higher 
accretion rates, the material seems to be dragged from more dis¬ 
tant places or, conversely, more distant parts reach the required 
temperature and density to produce emission lines (SA12). This 
is consistent with accretion models for the Ha line (Lima et al. 
2010). We do not hnd any signiheant correlation between the 
blueshifted and redshifted BC Gaussian components, which is a 
sign that they are independent and not caused by the same phys¬ 
ical process (for instance, an increase in the accretion rate trig¬ 
gering an accretion-powered wind). They seem only linked by 
temporal variability, as expected for bulk motion of material ro¬ 
tating as it falls onto the star. 

Finally, we searched for periodic signals in the BC. The cen¬ 
tering errors for a broad Gaussian are much larger than for the 
NC, so we also checked the location of the peak of the BC emis¬ 
sion (after subtracting the NC), which is less affected by the 
choice of ht. We only consider those data points where the am¬ 
plitude of the BC is >0.3 (over the normalised and subtracted 
continuum, well above the noise level). The results of the GLSP 
for the BC are listed in Table IA.2I For the BC Gaussian com¬ 
ponents, the periodicity analysis is not conclusive, mostly due 
to the small number of data points. The BC peak shows multi¬ 
ple low-signihcance periods, in agreement with the periodic or 
quasi-periodic signatures in the 7-8d range observed in the NC 
of the lines. Examining the shape of the BC line (see Appendix 
0, we also find a trend with the phase of the NC/photospheric 
absorption lines, especially for the cases with strongly bluesifted 
BC. This suggest that there is a common mechanism causing the 
velocity modulations of all emission lines and that the NC and 
BC are physically related. 


4. Dynamical and physical models for the spot, the 
star, and the accretion column 

In this section, we construct simple models to explain the RV 
signatures of the emission lines in EX Lupi, and explore the 
possibilities of a stratihed 3-D hot structure to reproduce the dy¬ 
namics and the strengths of emission lines. We assume that the 
footprint of the accretion structure behaves like a hot spot. A hot 
spot on the stellar surface would cause a clear modulation of the 
photospheric absorption lines as shown in Eigure 7j and could 
mimick the RV signature caused by a companion. The main chal¬ 
lenge for the spot scenario is that different lines have different 
RV modulations, which is inconsistent with plain flat hot spot 
models. Nevertheless, accretion can alter the vertical tempera¬ 
ture structure of the photosphere, suppressing normal absorption 
line formation in the accreting region. The extent to which the 
absorption is suppressed will depend on the strength and depth of 
the formation of the lines concerned (one example being the Li I 
6708 A line, the only photospheric line that still presents a small 
absorption component during outburst; SA12). The diversity of 
the amplitudes and the zero-point offsets between lines is a sign 
of a more complex 3-dimensional (3-D) accretion structure, ex- 


Table 1. Amplitudes and offsets of a sinusoidal fit to the emis¬ 
sion line signatures. 


Line 

A 

(km/s) 

w 

(km/s) 

<pQ 

(deg) 

He II4686 

3.06±0.28 

0.66±0.19 

11.8±0.8 

He I 4713 

1.74±0.26 

4.55±0.19 

3.1±1.4 

He I 5875 

1.03±0.19 

4.46±0.14 

-0.3±1.8 

He I 6678 

0.81±0.12 

3.28±0.09 

-2.7±1.5 

Fell 4351 

1.06±0.16 

0.70±0.11 

9.8±1.3 

Fe II4549 

1.10±0.18 

2.85±0.12 

8.5±1.4 

Fell 5316 

0.99±0.16 

1.86±0.11 

8.0±1.4 

Fe II 5362 

0.75±0.15 

1.26±0.11 

-0.19±1.9 

FeI4215 

1.05±0.16 

3.14±0.10 

4.8±1.2 

Fe I 4533 

1.58±0.33 

3.00±0.24 

0.7±2.0 

Fe I 5270 

0.86±0.14 

0.12±0.09 

6.2±1.4 

Fe I 5370 

1.50±0.24 

0.72±0.19 

2.9±1.7 

FeI5371 

1.15±0.18 

-0.06±0.12 

4.4±1.3 

Mg I 5167* 

0.83±0.17 

3.44±0.12 

4.8±1.8 

Mg I 5172 

0.71±0.13 

0.21±0.09 

4.8±1.3 

Mg I 5183 

1.06±0.14 

0.41±0.10 

3.9±1.2 

Ca II 8498^ 

0.31±0.07 

0.33±0.05 

-2.6±2.3 

Ca II 8662^ 

0.14±0.10 

1.46±0.07 

1.7±6.1 


Notes. Fitting the dynamical signatures (Figure with a sinusoidal 
modulation (RV=Vo-l-V,-l-Asin[2;r(())-l-(/io)]. ”A” (amplitude) refers to 
the strength of the sinusoidal modulation, ’’V^” (zero-point velocity) is 
the general offset of the line, after taking into account the zero-point ve¬ 
locity of the star (Vo=-0.52 km/s; K14), The phase is ”0” (for a 7.417d 
period), and ’Vo”, the phase offset. The values here correspond to the 
’’best fit”, but some of the parameters, such as the amplitude, have a 
relatively large spread (0.5-1.0 km/s) comparing data taken at different 
dates and similar phases. * The Mg I 5167A line may be affected by a 
nearby Fe I line. ^ The modulation of the Ca II lines is not significant, 
as explained in the text. The fits for the Ca II lines are listed to show the 
evident differences (especially, regarding the amplitude) with the rest of 
lines that have periodic modulations consistent with those of He II and 
the photospheric absorption lines. 


tended in the vertical (radial) direction and having a temperature 
and density gradient within the flow and/or accross the accretion 
column. This is consistent with the observations during outburst 
(SA12), and with evidence found in other T Tauri stars (Dupree 
et al. 2012; Petrov et al. 2014). 

4.1. The dynamics of narrow and broad components and 
their RV signatures 

The location of the line-emitting region can be inferred by mod¬ 
eling the observed dynamical signatures and radial velocities 
(Figure]^. If the NC lines are produced in a localised accretion 
spot structure, stellar rotation would affect their radial velocities 
imprinting a sinusoidal modulation. A similar scenario has been 
proposed for UV lines, for which some stars are also consistent 
with accretion-related spots getting in and out of view as the 
star rotates, and showing different lines-of-sight depending on 
whether they form in the accretion column or at the associated 
spot (Romanova et al. 2004; Ardila et al. 2013). Simulations pre¬ 
dict that the coverage of the spot decreases as the density of the 
accretion column increases and at small misaligment angles be¬ 
tween the stellar rotation and the magnetic moment (Romanova 
et al. 2004). Given that the accretion structures in EX Lupi are 
expected to be very dense to produce the observed emission lines 
(see Section [A^ , accretion models would favour very localised 
spot(s). The NC emission lines present a sinusoidal modulation 
with the same periodicity as observed for the absorption lines 
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Fig. 7. A scheme of how the photospheric absorption lines are affected when a hot spot moves in and out of our line of sight when 
the star rotates. In A we show the different parts of the stellar photosphere that produce the photospheric absorption lines. B shows 
how the line is affected when a hot spot is included. When the absorption deficit associated to the spot appears blueshifted, the 
centre of the photospheric lines seems to move towards the red. When the spot appears redshifted, the centre of the absorption line 
appears shifted towards the blue. Small emission lines are visible within some of the photospheric absorption lines, tracking the local 
rotational velocity of the part of the stellar surface where the line originates. Non-photospheric lines or lines that are particularly 
strong on the spot could appear simply as emission lines over the continuum, without any underlying absorption. 


(Figure]^, but variable amplitude and zero-point shift (Table[T]l. 
For lines with lower excitation potentials, the significance of the 
periods is lower and the signature is quasi-periodic rather than 
purely periodic. This can happen if they are produced in more 
extended and potentially messy structures (Figure]^. 


The amplitude is rather constant within a given species (e.g. 
He I, Fe I, Fe II), although there are significant shifts in the zero- 
point velocity, which is always redshifted. The maximum shift 
is observed for the He I lines, which is similar to the results of 
Beristain et al. (2001) on a large sample of CTTS, although the 
fact that this line could be in part affected by a redshifted absorp¬ 
tion (due for instance, to wind; Section |3.2| i could also led to a 
higher redshift. Differences in zero-point velocity can also result 
from differences in the optical thickness (and thus depth) of the 
lines (Section [3.2.1 1 . The differences in profile/velocity between 
high- and low-temperature lines would be associated to different 
parts of the accretion structure. Detailed models (Romanova et 
al. 2004, 2011; Ardila et al. 2013) have suggested various possi¬ 
bilities, ranging from vertical stratification (e.g., in an accretion 
structure with an aspect ratio where photons can easily escape 
through the walls), to horizontal stratification accross the ac¬ 
cretion structure (where the low-density, cooler material at the 
edges of the accretion column would also show a slower-than- 
free-fall motion). A further possibility, the presence of several 
accretion columns with various densities around the stellar sur¬ 
face, appears unlikely for EX Lupi, given that all lines show a 
very similar modulation (period and phase). An exception could 
be the Ca II IR NC lines, which do not show modulation and 
could be produced in a more uniform way over the stellar sur¬ 
face. 


We can use the observed sinusoidal modulation to investi¬ 
gate the typical location of the emitting structure (Romanova 
et al. 2004). The post-shock region at the bottom of the accre¬ 
tion column would be similar to a plage as seen in active stars 
(Ardila et al. 2013; Dumusque et al. 2014), although in this case, 
our ’’plage” would be 3-D and stratified in density, velocity, and 


temperature. If a plage (or accretion structure) is at a low lati¬ 
tude, it would only be visible during about half of the rotational 
period, producing an incomplete sinusoidal modulation (from 
blueshifted to redshifted; the part from redshifted to blueshifted 
would not be visible). We observe a full modulation, which re¬ 
quires the plage-like structure to be visible during the whole ro¬ 
tational period. If we observe the star under an inclination angle 
/, this happens for latitudes higher than i (or, in spherical coor¬ 
dinates, 6i9Q-i). The observed modulation would be caused by 
rotation and the 7.4Id period would correspond to the rotational 
period of the star. 

For an emission line produced in a small spot located on the 
surface of the star at spherical coordinate 0s, the rotation of the 
star will change the central velocity of the line following; 


RViine - Vq sin/ sin6(iSin(^ = Vq -Hcsin/ sin0jsin0(l) 

Here, Vq is any zero-point velocity that the system or the emit¬ 
ting zone may have (including the stellar velocity, Vo, and any 
systemic velocity in the accreting structure, such as a system¬ 
atic redshift in an infalling component, in Table [T]), P is the 
rotational period of the star, corresponds to the phase angle 
(from 0 to In), and nsin/ is the projected rotational velocity of 
the star. If the spot does not change its latitude during the time 
we observe it, we obtain a sinusoidal modulation with a radial 
velocity amplitude that depends on the viewing angle / and the 
spot location Og. A purely polar spot (Os-O', such as the one pro¬ 
posed by Dupree et al. 2012 for TW Hya) suffers no modulation. 
For an equatorial spot, the radial velocity amplitude is maxi¬ 
mal, although it will not be visible at all times. The maximum 
radial velocity that can be observed in the emission lines pro¬ 
duced at the surface of the star is thus usin/. Lines originating 
at some distance (d) over the surface of the star could have ra¬ 
dial velocity amplitudes up to 2;7r(R*H-d)sin //P, assuming that the 
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Fig. 8. A cartoon of a star with an accretion column that is visible at all times due to the stellar inclination and whose emission lines 
suffer Doppler shift as the star rotates (the observer is located perpendicular to the plane of the paper). The blue (red) column shows 
the position at which the lines originated in it would appear blueshifted (redshifted) to an observer looking from out of the page. 
The hashed areas mark the shock region. We suggest that the observed narrow emission lines are produced in the post-shock area, 
which is stratified in density and temperature. Lines produced in the shock region (He II, requiring very hot and dense conditions) 
originate in a nearly flat structure, so their periodic signatures are much stronger and less messy than other lines produced in the 
extended, infalling post-shock region. The BC is produced in the accretion column before the pre-shock region, showing higher 
velocity offsets as they rotate with/around the star. Not to scale; note that the shock region is at a very small distance over the 
photosphere, compared to the stellar radius. 


line-emitting structure rotates as a solid with the statj^ Similarly, 
if the extended structure where the lines are produced is dis¬ 
placed off the radial direction (for instance, trailing the star, or 
extended along the stellar surface), we would expect to observe 
some phase shift between lines originated at various altitudes. 
Although signihcant differences in zero-point velocity and phase 
are found comparing various lines (Table[^, optical depth differ¬ 
ences between lines and unresolved absorption (e.g. wind com¬ 
ponents) could also produce changes in both measured quanti¬ 
ties that cannot be properly estimated without detailed radiative 
transfer models. 

Assuming that the period is 7.417d and t;sin/=4.4 km/s, the 
system inclination angle i would be 23.8-25.5 degrees (for radius 
1.5-1.6Ro). For the He II line, we observe a mean amplitude of 
about 3 km/s, corresponding to 6s ~43 degrees (latitude ~47 
degrees). For an inclination of 23.8-25.5 degrees, this is compat¬ 
ible with a spot that is visible at all times, in agreement with the 
observed modulation. 

The information in Table [T] can be used to infer the rela¬ 
tive location of the various species along the emitting region, 
although for species with a large range in excitation conditions 
and critical densities (e.g. Fe I and Fe II), line strength and opac¬ 
ity differences can make this stratification not so evident. We 
explore two options for stratification; along the extended spot 

^ This may not hold as we move away from the object. Note that 
from the of the lines, if infall is assumed, they would be produced 
very close to the stellar surface, so d«R,. 


vs across the spot (Ardila et al. 2013). In the hrst case, material 
closer to the shock region would be expected to be hotter and 
denser compared to material closer to the stellar photosphere, 
and the relative velocities would depend on a mixture of infall 
and magnetic/gas pressure. In the second, the edges of the post¬ 
shock structure would be cooler and less dense, and simulations 
predict it to show slower infall velocities compared to the hot 
and dense column centre (Romanova et al. 2004). This situation 
is analog to the observations of UV lines in CTTS (Ardila et al. 
2013). 

The He II line has by far the strongest requirements regard¬ 
ing high temperature and density, which would require forma¬ 
tion relatively close to the shock and/or deep within the post¬ 
shock structure. In both cases, this spatial restriction would 
favour a stronger rotational modulation, compared to other lines, 
which is also observed. Depending on the velocity structure in 
the shock/post-shock region, its null zero-point velocity shift 
could suggest excitation very close to the shock and with a strong 
magnetic/gas pressure support accounting for the lowest infall 
velocity. In case of stratihcation accross the column. He II would 
be expected to form in the innermost part, although in this case 
its expected infall velocity would be higher than for lines formed 
in the periphery (Romanova et al. 2004). The difference in phase 
could be caused by location higher on the post-shock structure 
(compared to the stellar photosphere) or by a bow-shaped spot 
where the density varies along the stellar longitude (Romanova 
et al. 2004). 
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The He I and metallic lines would be produced in cooler, 
less dense parts of the post-shock region. Their small redshifted 
velocities suggest that the material, despite being basically at 
rest (due to magnetic/gas pressure), is still accreting or infalling. 
Taking into account the phase differences, after He II, most of 
the Fe II emission would be produced closer to the shock area or 
to the denser part on a bow-shaped spot (having the next large 
phase difference and lower infall velocity). He I lines also re¬ 
quire hot and dense material. The nearly zero phase difference 
observed for He I line and the relatively large infall velocities 
(if blueshifted absorption due to wind is negligible) would in 
fact suggest an origin very close to the stellar surface, although 
the higher He I abundance (compared to metals) can also result 
in emission over larger areas/volumes, including more extended 
parts over the stellar surface and eventually, winds. Fe I and Mg I 
lines would be mainly emitted in an area with lower temperature 
conditions, comparable to photospheric temperature/conditions 
of G/F stars. For an infall velocity of 4.5 km/s (typical for the He 
I lines), the material should be falling from 1.4x10 (~160 

km) over the stellar surface. Although gas/magnetic pressure 
and turbulence counteract the infall, models suggest that fully 
ionised post-shock columns have sizes of the order of 10^-10^ 
km (Dupree et al. 2012; Ardila et al. 2013), compatible with our 
findings of a relatively compact region compared to the stellar 
radius. Figure[^shows a scheme of the proposed configuration. 

We can also use Equation [T] to investigate the effects of such 
emission lines on the radial velocity of the absorption lines. We 
can model a photospheric absorption line using a Voigt pro¬ 
file that reproduces a typical, unblended photospheric absorp¬ 
tion line. We then add a small emission component with a cen¬ 
tral velocity that has a sinusoidal modulation such as shown in 
Equation[2 A continuum veiling contribution can also be added. 
We estimate the apparent radial velocity variation of the result¬ 
ing absorption line by cross-correlation. The modification of the 
stellar photosphere by the proximity of an accretion hot spot can 
be very complex, depending on plasma parameters, and can also 
contribute to powering a stellar wind (Orlando et al. 2010,2013). 
Therefore, this analysis does not pretend to derive strong con¬ 
straints on the properties of the line-dependent veiling, since it 
involves too many free (or poorly constrained) parameters, but to 
demonstrate that a simple accretion-plus-rotation model can ex¬ 
plain simultaneously the observed RV signatures in the emission 
and the absorption lines. 

We first considered the effect of the narrow emission lines, 
finding that although the presence of a small NC within the line 
can produce radial velocity offsets as large as observed, they also 
produce a very strong distortion of the line profile, which is not 
observed (some lines have narrow emission-line cores, but they 
were excluded from the radial velocity analysis in K14). We 
then consider the effect of a broad component in emission on 
the absorption line. Although broad components with emission 
over the continuum are only seen in certain lines (e.g. Ca IIIR 
triplet, strong Fe II lines such as the 5018 and 4923 A lines), we 
also see clear evidence of line-dependent veiling in EX Lupi, as 
proposed by Dodin & Lamzin (2012) for other CTTS. A single 
veiling factor (or a slowly changing one, such as resulting from 
a black-body approximation to the veiling) cannot reproduce all 
the lines, even if we restrict ourselves to a very small portion of 
the spectrum (see more details in Section [42| ). This is a sign that 
many of the absorption lines, if not all, have a line-dependent 
veiling, or extra emission that fills part of the absorption line 
and that is at least as broad as the absorption line itself (to avoid 
distorting the line profile). As the accretion rate increases, the 
line-dependent veiling would increase, filling in the whole ab¬ 


sorption line and eventually showing as a broad emission profile 
for very high accretion rates, as observed in outburst (SA12). 

Figure shows that, for reasonable values for the BC width 
and amplitude, it is possible to reproduce the observed radial ve¬ 
locities while avoiding a significant distortion of the line profiles 
and symmetry. Our simple model assumes that the rotationally- 
modulated BC in emission has a Voigt profile (with various 
Gaussians and Lorentzian contributions, y). A Gaussian core 
about 3-4 times broader than a typical photospheric absorption 
line avoids introducing large distortions in the line profile. We 
also consider a BC strength that results in veiling factors for 
the line between 0.25-0.65, and that it has a peak-to-peak ve¬ 
locity amplitude for in the range of 10-30 km/s and a few km/s 
net redshift (as observed; Table[T]). A velocity amplitude of ~20 
km/s would correspond to a distance of about ~5.7 stellar radii 
(assuming solid-body rotation and Os ~43 degrees), compatible 
with material near the edge of the stellar magnetosphere. We 
also include a ’’classical” continuum veiling as a constant emis¬ 
sion factor that dilutes both the photospheric absorption line and 
the wings of the emission line. Blending of nearby lines (which 
often occurs in the blue part of the spectrum) could also result 
in a relatively flat extra emission over the real continuum. The 
broad line causes a rotationally modulated ’’tilting continuum” 
that changes the apparent radial velocity of the line. With the 
strongest change being in the broad line wings that always fall 
out of the absorption line, the line bisector distortion is minimal 
and comparable to the distortion observed in K14. Checking the 
two bisector parameters explored by K14 (the bisector velocity 
BVS and bisector curvature BC, see K14 for details), the varia¬ 
tions observed in our model are comparable to the observations 
of the photospheric absorption lines. 

This simple model can reproduce a radial velocity modula¬ 
tion very similar to the one observed. A small net redshift and/or 
a non-Gaussian profile (adding Lorentzian wings) can produce 
the observed red-blue asymmetry in the RV curve (interpreted 
as eccentricity in the companion scenario). Further asymmetries 
in the observed curve could be caused by the emitting region not 
being entirely symmetric|^(due to occultation effects, or to vari¬ 
ations in the observed projected velocity for an extended spot as 
the star rotates), or by having a small phase offset (which would 
happen if the accretion structure is trailing the star or if the spot 
is extended along the stellar surface). As long as the period of 
the line-dependent veiling is maintained, small variations in the 
strength of the lines would not destroy the global modulation. 

We thus conclude that the most complete explanation of the 
observed RV modulation of the photospheric absorption and the 
emission lines is offered by assuming a modulation induced by 
line-dependent veiling throughout the absorption-line spectrum. 
It also shows that the phenomenon of line-dependent veiling by 
broad lines may be one of the worst-case scenarios to detect the 
presence of companions in accreting stars, given how robust the 
signal is to small variations in the accretion rate, and how weak 
the induced line asymmetries can be if the BC is significantly 
broader than the photospheric absorption lines, being diluted 
with emission from nearby lines and/or extra ’’classical” con¬ 
tinuum veiling. If most photospheric lines are affected by line 
veiling, distinguishing this effect from any other periodic mod¬ 
ulation (such as the one induced by a companion) may be very 
hard. 

This model requires that the accretion-related structure is 
stable over several years. The signature does not seem affected 


Note that for instance, the Ca II IR BC are not entirely symmetric 
but deviate from pure Gaussian profiles. 
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by the strong 2008 outburst (although there are only 3 points 
taken before the outburst), suggesting that the accretion struc¬ 
tures infalling onto this star are extraordinarily well-hxed to the 
stellar structure/magnetic field. Even though the accretion rate 
varies in time and is probably very clumpy, it would be a case of 
remarkably stable accretion columns with very little instability 
(Kurosawa & Romanova 2013). Checking the long-term stability 
of the signal will be also a strong test for the proposed accretion 
column scenario. 


4.2. Physical conditions in the accretion column 

We also explore whether the observed mean spectrum (obtained 
by combining all the normalised observed spectra) can be repro¬ 
duced as the sum of a stellar template, a classical veiling contin¬ 
uum (constant over small wavelength ranges) and emission lines 
originated in a region with a given temperature (T) and electron 
density (n^). Dynamics and BC emission are not considered in 
this simple model. We do not attempt to model the H Balmer 
series, Na I D, Ca IIIR triplet or other broad lines (or BC) since 
they are very complex and self-absorbed. We also do not have 
high enough S/N to constrain line ratios and physical condi¬ 
tions for the BC (which was nevertheless possible during out¬ 
burst; SA12). Lines with different excitation potentials and criti¬ 
cal densities allow us to explore the structure of the line-emitting 
accretion column in quiescence (S12; Petrov et al. 2014). 

For the photospheric model we consider an effective tem¬ 
perature Teff—3150 K, logg=4.0, and t;sin/=4.4 km/s (Sipos 
et al. 2009), using the standard template from Coelho et al. 
(1995). The veiling is considered as a continuum with a certain 
value over a small wavelength range. The atmospheric model 
with the veiling continuum should reproduce the observed ab¬ 
sorption spectrum, but we hnd that a simple continuum veiling 
with a slow dependency on the wavelength cannot reproduce si¬ 
multaneously all the photospheric lines, even if we only con¬ 
sider a small wavelenght range (Figure [T0|i. This suggests a line- 
dependent emission component hlling in part many (if not all) 
the emission lines with different strengths. 

For the emission line spectrum, we consider only the metallic 
and narrow lines, including He I and He II, which are dominated 
by narrow components. We hrst estimate whether the line emis¬ 
sion is optically thin or thick, and obtain approximate values (or 
ranges) for the temperature and electron density in the emitting 
region. In broad lines, the high velocity gradients ensure that 
the lines are optically thin (Beristain et al. 1998; Petrov et al. 
2014). For EX Lupi, the NC the line ratios are nevertheless con¬ 
sistent with partially optically thick emission. We use the atomic 
data from the National Institute of Standards and Technology 
database (NIST; Kelleher et al. 1999, Ralchenko et al. 201(Q to 
estimate line ratios in the optically thin approximation, which are 
proportional to the atomic gf factors for permitted lines (Martin 
& Wiese 1999). For a transition with upper level k and lower 
level i, the line strength g,f,i can be estimated as: 


^ifik 


lire 


( 2 ) 


where g, and gk are the level multiplicities (2J-I-1), A^, is the 
transition probability, c is the speed of light, m^, and e are the 
mass and charge of the electron, and eo is the dielectric constant. 
Comparing line pairs with similar upper-level energy (Table 
see Table A. 3 for a complete list of unblended lines) shows that 


http:physics.nist.gov/PhysRefData/ASD/lines_form.html 


although some of the lines (e.g. He I) are optically thick and sat¬ 
urated, the metallic neutral and ionized lines (Mg I, Fe I, Cr I, 
Mn I, Fe II, Cr II) span a broad range where the stronger lines 
tend to be optically thick and the weaker ones are sometimes 
optically thin. We also observe significant differences compared 
with the case of DR Tau in Beristain et al. (1998), such as the 
suppression of the lines with lower transition probabilities (A,y) 
and lower excitation temperatures, which point to higher temper¬ 
atures and densities and/or stronger high-energy irradiation for 
EX Lupi than for DR Tau, despite its lower accretion rate. 

The density can be estimated from the accretion rate (SA12): 


M 

Afv^i 


(3) 


Here, n is the number density, M is the accretion rate, A is the 
area of the stellar surface of which a fraction f is covered by 
the accretion-related hot spot, v is the velocity of the accreted 
material as it arrives to the shock region, and fj. is the mean 
particle mass. If we consider the observed values for the accre¬ 
tion rate (M~10“'°-10“® Mo/yr) and the velocity at which the 
material arrives to the shock region (v= 100-200 km/s), we ar¬ 
rive to n=6-120/f xlO^ cm“^. We assume that the spot covers 
a small part (~1%) of the stellar surface (suggested by Grosso 
et al. 2010 based on X-ray observations, and expected in case 
of higher density; Romanova et al. 2004). At relevant tempera¬ 
tures all the metals would be ionised and the H and He would 
be mostly neutral, producing electron densities of n^ ~ 10^-10'° 
cm~^. If the material arrives to the shock region at high velocities 
(v=100-200 km/s) but is deccelerated onto the star at few km/s 
(Table in the post-shock region where NC are originated, we 
can expect the density to increase by up to two orders of mag¬ 
nitude, arriving to n^ ~10'°-10'^. These higher values are more 
in agreement with what we would expect from the nearly equal 
ratios observed for the Ca II IR triplet, among others. Extra ion¬ 
ization due to energetic radiation in the accretion shock or in the 
stellar chromosphere could also increase the electron density in 
the post-shock area. The presence of extra sources of ionization 
were evident during the outburst (SA12), and some extra con¬ 
tribution of ionising UV emission (accretion-related) and by the 
shock itself (Orlando et al. 2010) is also likely to happen in qui¬ 
escence. 

We also follow Hamann & Persson (1992) to evaluate the 
electron density based on the requirements to achieve optically 
thick emission lines in the spectrum (similar for chromospheric 
emission and an accretion column or post-shock region sce¬ 
nario). The condition for equal flux in a multiplet of emis¬ 
sion lines (such as the Ca II IR triplet, or the Mg I triplet) 
is satisfied when collisional decay dominates over the transi¬ 
tion rate, rieCki »Aj:,7t. This mean Tn^, > 10'^ cm“^ for Ca II 
8542/8662A. Grinin & Mitskevich (1988) suggest that values of 
T >1 for the 8498A line, t >10 for the 8542A line are enough 
for this. Considering the lower r values. Shine & Linsky (1974) 
estimate the minimum threshold density for detection of Ca II 
emission to be Nca// ~ 2.5xl0'^(V/50 km/s) cm“^. Assuming a 
solar composition and that most of the Ca is ionised, this results 
in N(h/+h//) ~ 2.5x 102'(V/50 km/s) cm^^ 

In EX Lupi, the BC of the Ca II IR triplet becomes unde¬ 
tectable when the accretion rate is at its minimum (M~ xl0“^° 
Mo/yr). This suggests that the Ca II BC is very close to its opti¬ 
cally thin limit, so that n,, ~10^^-10'^ cm“^ in the region where 
the BC originates. If we assume that the accretion rate is constant 
for the amount of matter that flows through the BC-region into 
the NC-region, M=nBcA fsc fvc Vvc. the relation 
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Fig. 9. Results of our toy model for the radial velocity signatures (upper panel), and bisector velocity shift (BVS) and curvature 
(BC; lower panels). It shows the kind of radial velocity modulations that can be induced on a photospheric absorption line by an 
emission component (broader than the absorption line) that shows a rotational modulation. ’P’ denotes the peak value of an emission 
component with a Voigt profile with Gaussian width cr and Lorentzian contribution y. ’Cont’ marks the classical veiling continuum 
added to the model. ’A’ denotes the radial velocity amplitude of the emission line (maximum RV offset for a sinusoidal modulation), 
and ’RS’ is the zero-point redshift of the line. The dashed line box shows the values of the bisector velocity as measured by K14, 
which are similar to those derived from the toy model. Although there is a slight BVS slope, the typical errrors at this level can 
make it very hard to distinguish. 


between the density depends on the velocity ratios of the mate¬ 
rial, and the areas covered by the BC and NC zones. If there are 
no different sources of extra ionization between the two regions, 
the relation can be transformed into a relation for the electron 
density: 


ne,NC - ne,BC 


vbc 

VnC 


fBc \ 
Inc I 


(4) 


The velocity widths are about 12-14 km/s for the NC and 50-100 
km/s for the BC, which means by itself an increase in density by 
a factor of 3.5-8.5. If in addition the accreting structure gets thin¬ 
ner/compressed as it moves into the shock/post-shock region, a 
further increase in the density is to be expected. This means that 
in the NC emitting region, n^. > 5xl0'^-10'^ cm“^. We will use 
this estimate as a starting point in the subsequent calculations. 
Such density estimates are in agreement (slightly higher) with 
the values derived by Petrov et al. (2014). 

Following SA12, we apply Saha’s law for various tempera¬ 
ture and electron density values, in order to constrain the tem¬ 


perature and density structure in the accretion columns (Mihalas 
1978). Saha’s law works under the assumption that the system 
is in local thermodynamic equilibrium (LTE). The LTE assump¬ 
tion is questionable, although acceptable in an approximate toy 
model. We use the NIST line emission calculator to estimate line 
intensities and ratios for a given T, n^. As in outburst, a single 
temperature/density pair is unable to reproduce all lines at the 
same time. The lack of significant emission from Ti I lines puts a 
very strong constrain to the typical temperatures in the accretion 
column. Eor the relevant densities, a minimum temperature of 
~7000 K is needed to make the Ti I lines undetectable compared 
to Ti II. Regarding the He I/II emission, a temperature between 
8000-10000 K can explain the strength of the 5875 and 6678 A 
lines, but the formation of more energetic lines (such as the He 
14920A and especially the He II4686 A lines) require tempera¬ 
tures of the order of 20000-23000 K, which would render other 
elements such as Fe, Cr, and Mg fully ionised. 

For our final toy emission-line model we thus consider a con¬ 
stant electron density n(.=10'^ cm“^ and three different temper- 
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Fig. 10. An example of the attempt to reproduce the median EX Lupi spectrum (black line, with standard deviation per pixel marked 
in grey) by using a rotationally-broadened standard veiled by a continuum emission (magenta). The difference between the observed 
and veiled spectra is shown in pink. A constant veiling cannot reproduce the depths of all lines, even if we only consider a very 
small wavelength range. The differences are too large to be explained by the standard deviation or time variation observed between 
the different spectra. 


atures; 8000, 10000, and 23000 K (0.69, 0.86, and 2.0 eV, re¬ 
spectively). Note that in the real situation, both the electron den¬ 
sity and the temperature are expected to vary within the accre¬ 
tion column. Since most of the lines are optically thick, the line 
intensity is simply scaled to match the observations. The final 
spectrum is obtained by summing the photospheric component, 
a continuum veiling, and the lines. Figure[TT|shows the results of 
our simple model in several wavelength ranges. Despite its sim¬ 
plicity, our toy model confirms that a temperature stratification 
is needed to fit the observed emission. A range of temperatures 
and electron densities is also consistent with our naive picture 
in Figure with the most energetic He I and He II lines be¬ 
ing dominated by the high temperature component, closer to the 
shock region, and Fe II and lower-excitation He I lines being 
dominated by higher temperature emission than Mg I/Fe I lines. 
Further exploration would require a proper line radiative transfer 
modeling with detailed physical conditions along the accretion 
structure, which is beyond the scope of this paper. 


5. Discussion: The star, the spot, and the accretion 
coiumn 

5.1. A critical revision of the origin of the line and RV 
modulations 

The velocity modulations observed in EX Lupi emission lines 
were first described in SA12, concentrating on the most con¬ 
spicuous modulations (tens of km/s) observed during the 2008 
outburst. For the RV variations in the photospheric absorption 
lines, K14 proposed two possible origins; RV signatures induced 
by a BD companion (with msin /=14.7Mj„^ and semimajor axis 
0.063 AU), vs rotational modulation induced by possibly two hot 
spots on the stellar surface. Both the model of the BD companion 


and the model of ’’classical” hot spot(s) had problems to produce 
a complete explanation of the observations. Here we propose a 
3-D accretion column as the origin of the RV modulation of both 
emission and absorption lines. An accretion hot spot or accretion 
column is necessary to produce the temperature inversion that 
gives rise to the emission lines. Therefore, if accretion alone is 
able to explain all observations, it would offer the most simple 
explanation, without the need to add further elements. 

The main caveat of the companion model is that a 
companion-only scenario cannot explain why the emission and 
absorption line velocities are off-phase. A tentative solution 
could be sought if the accretion would mostly land on the 
BD companion or in a structure connecting both objects (e.g. 
Artymowicz & Lubow 1996). In that case, we would detect the 
photoshere of the bright primary (a MO star, the BD photosphere 
would not be detected due to contrast) and the hot spot on the 
secondary or on the connecting structure. This would produce an 
anticorrelation between the RV of absorption and emission lines, 
similar to what has been observed in X-ray studies from accre¬ 
tion shocks in interacting binaries (Donatti et al. 2011, Argiroffi 
et al. 2012). Since the masses of EX Lupi and the proposed com¬ 
panion are very different, the amplitudes of the RV signatures 
for both would be very different. Nevertheless, this scenario is 
inconsistent with the observations. Considering the line with the 
strongest RV signature (He II 4686A line), there is only a very 
mild anticorrelation: a Spearman rank test gives a correlation 
coefficient r=-0.20 and false-alarm probability p=0.20, which is 
non-significant. Figure reveals that the emission line modu¬ 
lation is not simply the (scaled) stellar RV with a changed sign 
(which would be strongly asymmetric), but a redshifted symmet¬ 
ric sinusoidal modulation (Section |4.1| i. The timing of the RV 
cycle as observed in the absorption and emission lines is also in 
contradiction with the expected signatures in case of accretion 
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^ Observed data 

- ■ Total model spectrum 

— Standard, T,;j,^=3750 K, logg=4.0. 


— TiI/II0.69eV, n, = lel3cm 

— Hel/II 0.69eV, n,, = lel3cm 

— FeI/FeII0.69eV, n^ = lel3cm 

— Fel/Fell 0.86eV, n,, = lel3cm 

— FeI/FeII2.0eV, n^ = lel3cm •’ 
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Fig. 11. Toy model for the narrow emission lines, obtained by considering an electron density 10^^ cm~^ and three different tem¬ 
peratures along the accretion column (8000, 10000, 23000 K). The rotationally-broadened standard star spectrum is shown in grey, 
and the individual lines are shown in various colors. The final model, including a constant veiling, is obtained by summing all three 
components. 
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Table 2. Line ratios in the optically thin limit. 


Lines 

Observed 

Ratio 

Opt.Thin 

Limit 

Comments 

Hel 6678/5875 

0.41 

9e-5 

Very opt thick 

Mnl 4034/4033 

0.76 

0.64 

Nearly opt. thin 

Mnl 4034/4030 

0.72 

0.47 

Rather opt. thick 

Mgl 5173/5184 

0.92 

0.59 

Rather opt. thick 

Till 4313/4300 

0.53 

0.40 

Nearly opt. thin 

Till 4307/4302 

1.50 

0.76 

Optically thick 

Till 4054/4028 

0.63 

0.63 

Opt. thin 

Till 4450/4443 

1.0 

0.18 

Opt. thick 

Till 4590/4572 

0.78 

0.15 

Rather opt. thick 

Sill 6371/6347 

0.63 

0.59 

Nearly opt. thin 

All 3944/3961 

1.2 

0.50 

Opt. thick 

CrI 4290/4275 

0.80 

0.74 

Nearly opt. thin 

CrI 4290/4254 

0.73 

0.57 

Nearly opt. thin 

CrII 4617/4588 

0.43 

0.22 

Rather opt. thick 

CrII 4617/4559 

0.38 

0.23 

Rather opt. thick 

CrII 4876/4848 

1.0 

0.47 

Very weak lines 

CrII 4876/4824 

0.75 

0.58 

Very weak lines 

OI 7775/7774 

0.88 

0.60 

Rather opt. thick 

OI 7775/7772 

0.65 

0.43 

Rather opt. thick 

Fel 3896/3900 

0.68 

0.73 

Opt. thin 

Fel 3920/3928 

0.66 

0.60 

Opt. thin 

Fel 3923/3930 

1.07 

0.70 

Opt. thick 

Fel 4144/4045 

0.96 

0.16 

Opt. thick 

Fel 4216/4376 

0.82 

0.47 

Rather opt. thick 

Fel 5329/5172 

0.80 

0.88 

Contaminated by Mg I 

Fel 5456/5406 

0.82 

0.57 

Rather opt. thick 

Fell 4576/4351 

0.30 

0.15 

Rather opt. thick 

Fell 4520/4549 

0.39 

0.13 

Rather opt. thick 

Fell 4555/4549 

0.68 

0.30 

Rather opt. thick 

Fell 4583/4576 

0.76 

0.72 

Opt. thin 

Fell 4923/5018 

1.16 

0.73 

Opt. thick 

Fell 5235/5198 

1.16 

0.70 

Opt. thick 

Fell 6238/6248 

0.50 

0.47 

Nearly opt. thin 


Notes. Observed and optically-thin limit line ratios from similar-energy 
(calculated using the atomic data from the NIST database). 


in a binary system (Artymowicz & Lubow 1996), since it would 
request a (variable) difference in phase of approximately 1/4 of 
cycle, while accretion is most likely to occur along the line con¬ 
necting both objects. The phase variation between the compan¬ 
ion and the accretion would also require the accreted material to 
move at unphysical, variable speeds, which change the relative 
position between the spot and the companion during the phase. 
A scenario with spots on both objects would similarly fail, as it 
would either dilute the signature of the RV of the emission lines 
or modulate the line thickness (by summing a redshifted and a 
blueshifted components that change along the orbital phases). It 
also fails to explain why the RV amplitudes are often larger for 
the emission lines than for the photospheric absorption lines, and 
why the shapes of both curves are so different (asymmetric for 
absorption lines, symmetric for emission lines). 

The main caveat of the hot spot scenario in K14 is that it 
would require the 7.4Id period to be the rotational period of 
the star, which would be hard to reconcile with the v sin / and 
the constraints on inclination from disk models and CO obser¬ 
vations. We nevertheless hnd that usin/=4.4+2.0 km/s (Sipos 
et al. 2009) reproduces better the observations than the upper 
limit value of 2-3 km/s given in K14. With this higher rotational 
value, a rotational period of 7.4Id would require an inclination 
angle of 23-25 degrees (consistent with the disk models of Sipos 
et al. 2009). A low viewing angle was also predicted by Grosso 


et al. (2010) based on X-ray absorption. The observed veloc¬ 
ity amplitude of the NC is also in agreement with an origin in a 
high-latitude structure, in a star that is seen at a relatively low an¬ 
gle (Section [4/T] i. Moreover, all the NC velocities are redshifted, 
also in agreement with a relatively low viewing angle where the 
line-emitting region is always seen infalling (never blueshifted) 
as the star rotates. 

The problem for the hot spot hypothesis in K14 is that mod¬ 
els for a classical, continuum-emitting, flat hot spot would re¬ 
quire a coverage as large (or more) than the stellar surface to pro¬ 
duce the radial velocity signatures observed. Such a spot would 
also produce a luminosity modulation that is not observed in the 
photometry records. K14 assumed that the spot would suffer oc- 
cultations as the star rotates, which would increase the luminos¬ 
ity contrast, but a high-latitude spot visible at all times would 
not induce dramatic variations in the luminosity. In addition, 
if most or a large fraction of the energy released in the shock 
structure is emitted in lines (including many metallic UV lines) 
and not in the optical continuum, photometric rotational modu¬ 
lations may not be evident, while RV modulations would be sig- 
nihcantly enhanced with respect to the continuum-only scenario. 
This situation has been described in other T Tauri stars (usually, 
with strong accretion) and named ”veiling-by-lines” (Dodin & 
Lamzin 2012, 2013). 

The detailed line-by-line dynamical analysis breaks one of 
the main counter-arguments in K14 against the spot/column ex¬ 
planation, that the RV shifts observed in the emission lines had 
a smaller amplitude than the RV of the photospheric absorp¬ 
tion lines and that the signal periodicity was not significant. 
This happened because the relative offsets between lines and 
the various degrees of modulation (for instance, comparing He 
II with Ca II NC) dilute the global signature when many lines 
are combined. Once the lines are treated individually (includ¬ 
ing detailed individual fits), the modulation is strong and sig- 
nihcant. This is an extension of what we observed during out¬ 
burst (SA12), with different lines having various velocity off¬ 
sets depending on their excitation and physical conditions in 
the extended non-axisymmetric accretion structures. Petrov et al. 
(2014) also demonstrated that the stratification in velocities and 
physical conditions in the accretion columns result in different 
velocities for lines with various excitation potentials, although 
in this case the lines were seen in absorption and were not time- 
resolved. 

If the emission lines are formed in a hot spot near/on the 
stellar surface whose velocity is modulated by the rotation of the 
star, the largest extent of the spot is visible at zero velocity, so 
the line would be stronger when the spot is facing us. Despite 
the substantial variations from day-to-day and epoch-to-epoch 
in the line strength, there is a low-significance modulation of the 
amplitude/strength in several Fe II, Fe I, Ti II and Mg I lines. 
The strength of the lines is smaller at zero velocity, and maxi¬ 
mum when both the redshifted and blueshifted velocity offsets 
are maximal. Although this is the opposite situation than above 
mentioned, it is similar to the case of rotationally modulated X- 
ray emission (e.g. Flaccomio et al. 2005; Argiroffi et al. 2011, 
2012). X-ray emission is produced at the bottom of a 3-D ac¬ 
cretion shock, which results in the accretion column obscuring 
most of it at zero velocity, while most of the emission should 
escape through the sides of the accretion column at maximum 
redshift/blueshift. This effect has been observed for EX Lupi 
(Grosso et al. 2010). The same explanation could apply to the 
NC lines, if formed within a vertically extended structure where 
the emission can escape better along the ’’sides” of the accre¬ 
tion structure, instead of from the top. The Ca IIIR lines lack 
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these occultation effects and tend to be stronger when the struc¬ 
ture is facing the observer, in agreement with their origin in a 
less embedded, more delocalised structure with lower density 
and temperature. 

The rotational modulation gets increasingly smaller as we 
go to lines that are produced at lower-density, lower-temperature 
conditions, such as the Ca II IR line (Figure |^. Since lower- 
excitation lines require less restrictive conditions for their for¬ 
mation that can be met over a larger area (more similar to the 
conditions at the photosphere of the star), the RV signature of 
these NC lines could be diluted and/or appear rather as a quasi- 
periodic modulation. Differences in the rotation in various parts 
of the structure are also expected from models, and could lead 
to quasi-periodic modulations as observed here for the low- 
excitation lines (Romanova et al. 2004). Deviations from pure 
solid-body rotation in an extended structure can also dilute the 
modulatiorj^ The lower temperature, lower density areas will 
be larger, diluting the signature with a broader range of veloci¬ 
ties. Although the periodicity of the BC of the Ca II lines (the 
only BC detected in most of the spectra) is not significant, it is 
also consistent with modulations in the 7-8d range. If the ma¬ 
terial is fed to the star in a non-axisymmetric way, as proposed 
in outburst (SA12) and conhrmed by CO observations (Goto et 
al. 2011; Banzatti et al. 2015), it could also help to keep accre¬ 
tion confined to a relatively small area on/near the stellar sur¬ 
face, compared to cases where the material approaches the star 
from all directions. An accretion column/structure that trails the 
star (or material that is both accreting and rotating) could induce 
phase shifts in the rotational modulation, which is also observed 
(Table[T). 

5.2. Origin and impiications of stabie accretion coiumns in 
EX Lupi 

Accretion in EX Lupi appears to follow a combination of the 
standard accretion column scenario as an origin for the emis¬ 
sion lines in CTTS (Muzerolle et al. 1998, 2001), and the post¬ 
shock scenario (Dupree et al. 2012). We hnd evidence of both 
broad components consistent with an origin before the shock re¬ 
gion in an extended accretion column, and narrow components 
produced in a post-shock, slightly infalling structure extending 
down to the stellar photosphere. The line widths observed in the 
NC would suggest a less turbulent environment than proposed 
by Dupree et al. (2012) for TW Hya, being close to the expected 
thermal broadening for the relevant temperatures. The velocity 
differences observed between BC and NC are also consistent 
with the differences between a pre-shock and post-shock region 
(Ardila et al. 2013). The BC lines require structures with a large 
velocity span, which could be due to spatially extended rotat¬ 
ing/infalling packages of gas, although line broadening contri¬ 
butions from turbulence cannot be excluded. The lack of corre¬ 
lation between Ha/Ca II IR and He I lines is in agreement with 
the time delay observed by Dupree et al. (2012), although we do 
not have enough time-resolved coverage to follow the accreting 
packages on a short (hours) timescale. 

Cooling in the post-shock region would be dominated by line 
emission (including the metallic and He lines observed). Given 
the observed line velocities and stratihcation, the possibility that 
the narrow emission lines are produced in the post-shock re¬ 
gion is more likely than an origin in active regions of the star. 


® The region where the NC are produced is expected to be very small 
compared to the stellar radius, which means that the nearly-solid-body 
rotation would not need to be extended over a large vertical structure. 


which would tend to be more distributed over the stellar surface 
(Ardila et al. 2013; Brickhouse et al. 2010), even though many of 
the lines observed in emission are also common in active stars. 
A strong contribution from distributed chromospheric NC emis¬ 
sion would also tend to dilute the rotational signatures observed. 
Although this does not happen for most of the metallic lines, it 
could affect some of the typical activity indicators such as the Ca 
II NC. At this point, we cannot thus determine whether EX Lupi 
is more or less active than other similar young stars. 

It also remains unknown why most of the shock energy in 
EX Lupi is released in spectral lines. Stars with very rich emis¬ 
sion line spectra have typically high accretion rates (Hamann & 
Persson 1992; Lang et al. 2009; Sicilia-Aguilar et al. 2010). This 
is not the case for EX Lupi in quiescence, since most MO young 
stars are accreting at levels higher than ~10“'° Mo/yr (e.g. Lang 
et al. 2009; Sicilia-Aguilar et al. 2010). A plausible explanation 
would be a significant structural difference in the accretion chan¬ 
nels, since the emission of many of the observed lines (Ee II, He 
II) requires particularly hot and dense material. A similar ac¬ 
cretion rate being channeled through and landing on a more ex¬ 
tended part of the stellar surface may thus not be enough to pro¬ 
duce the rich forest of lines we observe in EX Lupi. The compar¬ 
ison of EX Lupi and DR Tau (from Beristain et al. 1998) shows 
that despite its lower accretion rate, the accretion structures in 
EX Lupi probably reach higher densities and temperatures. A 
lower temperature contrast between the stellar photosphere and 
the accretion structure would also tend to prevent the formation 
of the high energy lines. Interestingly, some of the line-rich ob¬ 
jects identihed by Lang et al. (2009) also show, like EX Lupi, low 
2MASS near-IR excesses and large Spitzer mid-IR fluxes. This 
could suggest a correlation between the strong line emission and 
the presence of inner holes or optically thin inner disks, although 
the sample of objects is too small to draw any statistically sig- 
nihcant conclusion. In addition, only ASASSN13db (Holoien et 
al. 2014), also an outbursting star, has a comparable spectrum. 
Whether the rich emission line spectra reflect a particular type 
of accretion structures in these EXors, or is a consequence of the 
accretion variability of the stars, is to be explored in the future. 

If line veiling (instead of continuum emission) is dominant 
for a signihcant number of stars, accretion rates are likely to 
be overestimated (Dodin & Lamzin 2012). Veiling is the ulti¬ 
mate calibration of accretion rate tracers, including U-band ex¬ 
cess (Gullbring et al. 1998) and emission lines (Natta et al. 2004; 
Lang et al. 2009; Alcala et al. 2014). This could have led to an 
uncertain estimate of the accretion rate of EX Lupi, especially 
in quiescence. In outburst, the accretion luminosity overwhelms 
the stellar luminosity, so the measured rate should be more accu¬ 
rate. The outburst accretion rate (few times 10“® Mo/yr; SA12; 
Juhasz et al. 2012) is high, but not extreme, for a low-mass star, 
which is also in agreement with a very low rate during the qui¬ 
escence phases. 

The observed velocities of the NC and BC can be compared 
with the different velocities associated to the star, which can be 
due to stellar rotation (assuming a structure would rotate as a 
solid body together with the star), infall (considered as free- 
falQ, and Keplerian rotation. Table shows a small summary 
of velocities and associated locations for these three cases (note 
that not all velocities are expected to be relevant at all locations). 
Taking into account the period of the star, the corotation radius is 
located at ~8.45 R», or ~0.063 AU, corresponding to a velocity 


^ Note that in the presence of a strong magnetic field near the stellar 
photosphere, the infall is probably caused by downward motions along 
the magnetic field structure that corotates with the star. 
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of ~92 km/s. This velocity is not too different from the maximum 
span of the peak of the Ca II BC (~ +100 km/s). As we move 
away from the star, Keplerian rotation (probably combined with 
infall) is probably dominant. Solid-body rotation with the star 
(even if the column is trailing the star) probably would break at 
relatively short distances, depending on the strength and extent 
of the magnetic field. Considering the observations of BC during 
outburst (SA12) and the CO data from Goto et al. (2011) and 
Banzatti et al. (2015), the non-axisymmetric, rotating/infalling 
structure of gaseous material may extend well into the few tenths 
of AU within the inner disk dusty hole. Feeding the star in a 
non-axisymmetric way may help to keep the rotational modula¬ 
tion visible over time (for instance, bringing material only onto 
certain areas over the stellar magnetosphere). 

An open question is whether the magnetic field structure of 
a star like EX Lupi can be stable over such a long period of time 
(~5 yr). Although rapid changes in the magnetic field configu¬ 
ration have been proposed for young stars (Brown et al. 2011), 
some observational evidence suggests that magnetic field rever¬ 
sals are not so common in young (10-50 Myr), rapidly rotat¬ 
ing, solar-type stars, at least in timescales of 3-7 years (Marsden 
et al. 2009). Rotational modulation is observed in some cases 
(Johns-Krull et al. 2013), but there are large differences in the 
behaviour of different low-mass stars, with some showing strong 
variability, while others have persistent fields over several years 
(Symington et al. 2005). As a young, low-mass star, EX Lupi 
should be completely convective, favoring a simple dipolar field 
(Gregory et al. 2014), which could lead to a simple accreting 
structure within a well-defined part of the stellar surface, as pro¬ 
posed here. Young stars with simple magnetic fields tend to have 
large disk truncation radii and spots at high latitudes (Johnstone 
et al. 2014), both of which would be in agreement with the loca¬ 
tion of the accretion structure and corotation radius for EX Lupi. 
On the other hand, close-in companions can have very complex 
field structures, leading to accretion distributed over large sur¬ 
faces on the stars (Donati et al. 2011). Although very different in 
mass, a hypothetical BD companion near EX Lupi could have a 
strong field (Morin et al., 2010) and lead to a complex total field, 
which would wash out the velocity signatures in the accretion- 
related lines. There are no constraints on the long-term magnetic 
field variability of EX Lupi (nor of large samples of CTTS), but 
in principle a magnetic field stability over several years would 
not be exceptional. Euture studies of the magnetic field of EX 
Lupi and its variability would be a very strong test of the valid¬ 
ity of our stable accretion hypothesis. 

6. Summary and conclusions 

We analyse the rich emission line spectra of EX Lupi to inves¬ 
tigate the origin of the emission lines and whether they could 
be responsible for the RV signatures previously observed. Our 
study shows how we can turn a well-known problem, namely the 
impossibility to obtain reliable radial velocities for young stars 
with emission lines to search for companions, into a useful tool 
to trace the structure and properties of the accretion columns. 
This tool is able to probe into regions over the stellar photo¬ 
sphere that are smaller than the stellar radii, which are inacces¬ 
sible by other means, including interferometry. 

- As previously observed, the EX Lupi spectra in quiescence 
display a very rich collection of metallic lines, most of which 
are narrow, although some strong lines have a NC+BC struc¬ 
ture. The BC becomes more prominent during periods of in¬ 
creased accretion. 


Table 3. Different velocities associated to the system. 


Velocity 

(km/s) 

^rotation 

(R.) 

^freefall 

(R*) 

P 

^kepler 

AU 

5 

0.14 

1.7e-4 

21.3 

10 

1.3 

7.0e-4 

5.3 

20 

3.5 

2.8e-3 

1.3 

50 

10.4 

1.8e-2 

0.2 (29R.) 

100 

21.7 (0.17 AU) 

7.5e-2 

0.053 (7.2R.) 

200 

44. (0.34 AU) 

0.39 

0.013 (1.8R.) 

300 

67. (0.51 AU) 

1.7 

5.9e-3 (0.79R.) 


Notes. Estimated for a rotational period 7.417d, M,=0.6Mq, R, = 1.6Ro. 
’’Rotation” refers to solid-body rotation (locked to the star). The dis¬ 
tances are measured from the stellar photosphere at R,. Note that for the 
rotation options, the velocities are corrected by the sini. Also note that 
at any given location, not all three motions are expected to be relevant, 
(e.g. solid-body rotation is unlikely very far from the star; Keplerian 
rotation is unlikely very close to the star; infall near the star is proba¬ 
bly dominated by the magnetic field and thus not pure free-fall), so this 
table is merely informative. 

- The spectra show line-dependent veiling. If a faint BC partly 
fills in many (or most) of the photospheric absorption lines, it 
can account for a substantial fraction of the accretion-related 
released energy. If most of the veiling in EX Lupi is emit¬ 
ted in lines and not as a hot continuum, there could be some 
uncertainties in the accretion rate estimates. Eor line-veiling, 
the total energy released would be smaller than in the contin¬ 
uum case. Depending on whether the stars used for the cal¬ 
ibration are dominated by continuum or line veiling, the ac¬ 
cretion rate could be lower. Although veiling-by-lines tends 
occur in strong accretors (Dodin & Lamzin 2012; Ingleby et 
al. 2013), EX Lupi demonstrates that objects with low accre¬ 
tion rates can also have significant line-dependent veiling. 

- Line-dependent veiling can explain the origin of the re¬ 
markable spectrum observed during the accretion outburst 
(SA12). The BC-dominated outburst spectrum, which hides 
all photospheric features under a whealth of broad emission 
lines, would correspond to the reaction of the line-dependent 
veiling to a higher accretion rate. During quiescence, the BC 
of most lines would remain undetected as it does not com¬ 
pletely fill in the photospheric lines. The same effect could be 
expected also in the outbursting star ASASSN13db (Holoien 
et al. 2014), which has a similarly rich outburst spectrum. 

- The NC emission is consistent with an origin in a hot post¬ 
shock region, deep in the accretion channels and close to the 
stellar photosphere. The dynamical analysis suggests that the 
accretion ’’spot” is a 3-D structure similar to a stellar chro¬ 
mosphere, where the temperature goes through a minimum 
at the top of the stellar photosphere and then rises afterwards, 
although on a localised place near the stellar surface. The 
profiles of the stellar absorption spectral lines can have nar¬ 
row emission cores if they are opaque enough to be formed 
in the region where the temperature is increasing outwards. 
The BC is consistent with being formed at larger distances 
from the star in rotating/infalling non-axisymmetric accre¬ 
tion structure(s), in agreement with the outburst observa¬ 
tions. In quiescence, the emitting structure appears closer to 
the stellar surface than in outburst. 

- There is a clear difference in velocity modulation between 
the NC and the BC. The BC shows velocities of several tens 
to hundreds km/s, while the material associated to the NC 
only shows infall at a few km/s rate. There are also clear 
differences between the various emitting species, consistent 
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with stratification in temperature and density within an ex¬ 
tended, 3-D accretion structure. The analysis of the dynam¬ 
ics of lines with different excitation potentials gives a self- 
consistent view of the accretion column structure and can 
be applied to further objects in order to trace spatial scales 
which are beyond direct spatial resolution. 

- The RV signatures observed in emission and absorption lines 
can be both explained with rotationally-modulated signa¬ 
tures of an accretion column that cools mostly by emission 
lines. We show that the line-dependent veiling emitted by the 
accretion column, velocity-shifted due to the stellar rotation, 
can produce a RV signature in the photospheric absorption 
lines that mimics the signature of a companion. The result¬ 
ing ’’tilting continuum” fills in the photospheric absorption 
lines, leading to an apparent RV shift. The velocity structure 
and phased RV curve in the emission lines can be used to dis¬ 
tinguish companion-related from accretion-related phenom¬ 
ena in young stars. In case of very stable accretion structures 
and line-dependent veiling, it can be nevertheless very hard 
to distinguish both scenarios. 

- The line velocity modulations observed for EX Lupi are 
suggestive of the stable accretion scenario proposed by 
Kurosawa & Romanova (2013). Our long-term monitoring is 
consistent with accretion occurring through remarkable sta¬ 
ble channels, despite the frequent variations in the accretion 
rate. The only difference between the quiescence and the out¬ 
burst phases is the amount of matter channelled onto the star. 
In order to maintain stable accretion channels over years, a 
very strong mechanism (maybe magnetically controlled or 
controlled by the disk) has to rule the flow of material onto 
the star. Given that the gaseous material within the inner 
disk hole also consist of infalling/rotating, non-axisymmetric 
flows of matter (Goto et al. 2011; Kospal et al. 2011; SA12; 
Banzatti et al. 2015), part of the flow regulation could happen 
within the inner disk. 

- The cause of the outbursting behaviour in EX Lupi remains 
unknown. Its disk is not substantially different from the 
disks around other stars. We cannot conclude that the star 
is anomalously active. The only remarkable fact is the stabil¬ 
ity of the accretion channels, which should be explored in the 
coming years. The potential relation of ’’overregulated” sta¬ 
ble accretion onto the star and the violent outbursts and vari¬ 
ability is an intriguing possibility that needs to be explored 
in the future (either as a cause or consequence of the out¬ 
bursts, or as star-controlled irregular accretion). A mismatch 
between the rate at which the disk feeds material onto the in¬ 
ner disk and the rate at which the star lets accretion proceed 
could be a possibility to explain the outbursts and the obser¬ 
vations of Banzatti et al. (2015) of material accumulated in 
the inner disk. Extending the detailed metallic line emission 
analysis to other (non-outbursting) stars may reveal whether 
this is a typical behaviour of some CTTS, or whether it may 
have a deeper relation to the accretion variations observed in 
EX Lupi and other EXors. 
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Table A.l. List of observations during EX Lupi quiescence analysed in this paper. 


ID 

Instrument 

Accretion 

BC shape 

2454309.615 

EEROS 

M 


2454310.657 

EEROS 

M 


2454311.687 

EEROS 

L 


2454874.838 

HARPS 

H 


2454875.850 

HARPS 

H 


2454876.808 

HARPS 

M 


2454877.791 

HARPS 

M 


2454877.906 

HARPS 

L 


2454891.800 

HARPS 

H 


2454891.882 

HARPS 

H 


2454892.774 

HARPS 

M 


2454894.760 

HARPS 

L 


2455053.675 

EEROS 

L 


2455058.493 

EEROS 

M 

B 

2455257.820 

EEROS 

L 


2455259.877 

EEROS 

L 


2455262.807 

EEROS 

L 


2455263.868 

EEROS 

L 


2455309.894 

EEROS 

L 


2455338.791 

EEROS 

L 


2455353.827 

EEROS 

L 


2455396.589 

EEROS 

H 

C 

2455396.637 

EEROS 

H 

C 

2455397.611 

EEROS 

M 


2455399.658 

EEROS 

M 

C 

2455400.536 

EEROS 

H 

B 

2455401.506 

EEROS 

H 

B 

2455402.584 

EEROS 

M 


2455403.612 

EEROS 

M 

C 

2455405.661 

EEROS 

H 

C 

2455407.620 

EEROS 

H 

C 

2455586.871 

EEROS 

H 

D 

2455587.853 

EEROS 

H 

C 

2455588.858 

EEROS 

H 


2455630.816 

EEROS 

H 

B 

2455637.895 

EEROS 

M 

B 

2455667.701 

EEROS 

H 


2455668.691 

EEROS 

H 

D 

2455669.737 

EEROS 

H 

C 

2455670.794 

EEROS 

H 

D 

2455671.698 

EEROS 

H 

C 

2455672.876 

EEROS 

M 


2455673.849 

EEROS 

M 

R 

2455674.824 

EEROS 

H 

B 

2456111.513 

EEROS 

L 


2456113.565 

EEROS 

L 


2456114.477 

EEROS 

L 

R 

2456115.477 

EEROS 

L 


2456116.524 

EEROS 

L 


2456117.562 

EEROS 

M 


2456118.670 

EEROS 

L 


2456119.494 

EEROS 

L 


2456120.508 

EEROS 

L 


2456121.579 

EEROS 

L 



Notes. The "Accretion” column correspond to the classification as high/medium/low accretion given in the text, based on the strength of the BC 
of the emission lines (see Appendix]^. The ”BC shape” field indicates the shape of the Ca IIIR broad component, which can be classified as 
redshifted (R), blueshifted (B), centred (C), or double-peaked (D; see Appendix]^. Only spectra with strong BC can be classified (see text). 


Appendix A: Observations and summary tabies 
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Notes. Except when otherwise stated, the signature correspond to the line NC. The field ’N’ indicates the number of lines used in the GLSP estimation, after removing low S/N and anomalous ones. 
NSP stands for ”no significant period”. Energy includes excitation potential (E*) and (for ions) ionization potential (E,), which are shown separately as a sum for clarity. Transition probabilities and 
line excitation potentials are not available for some of the lines. All atomic data are from the NIST database. 
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Table A.3. Non-blended, narrow emission lines observed in EX Lupi. Based 
on line identifications by SA12, and using the atomic constants fron the NIST 
database. Only lines dominated by a NC and not blended with other emission 
lines are listed here. 


'^ohs 

(A) 

EW 

(A) 

FWHM 

(A) 

Spec. 

^theo 

(A) 

^ki 

S~^ 

Ei - Ei 
(eV) 

" J/: 

3964.71 

-0.20 

0.37 

Hel 

3964.73 

6.95eH-06 

20.61 - 23.74 

0- 1 

4026.28 

-0.53 

0.51 

Hel 

4026.19 

3.22eH-05 

20.96 - 24.04 

2- 1 

4471.57 

-1.06 

0.49 

Hel 

4471.48 

2.46eH-07 

20.96 - 23.73 

2-3 

4713.24 

-0.20 

0.44 

Hel 

4713.15 

5.29eH-06 

20.96 - 23.59 

2- 1 

4921.99 

-0.17 

0.28 

Hel 

4921.93 

1.99eH-07 

21.22-23.74 

1 -2 

5015.69 

-0.43 

0.35 

Hel 

5015.68 

1.34eH-07 

20.62 - 23.09 

0- 1 

5875.74 

-1.67 

0.74 

Hel 

5875.62 

7.07eH-07 

20.96 - 23.07 

2-3 

6678.23 

-0.69 

0.59 

Hel 

6678.15 

6.68eH-03 

21.22 - 23.07 

1 -2 

7065.30 

-0.49 

0.79 

Hel 

7065.19 

1.55eH-07 

20.96 - 22.72 

2- 1 

4685.84 

-0.59 

0.66 

Hell 

4685.3 

1.47eH-07 

48.37-51.01 

5/2 - 5/2 

4287.89 

-0.11 

0.29 

Nal 

4287.84 

— 

2.10-4.99 

1/2 - 1/2 

4226.72 

-0.46 

0.22 

Cal 

4226.73 

2.18eH-08 

0.00 - 2.93 

0- 1 

6456.86 

-0.04 

0.34 

Call 

6456.87 

— 

— 

— 

8248.84 

-0.05 

0.38 

Call 

8248.80 

6.10eH-07 

7.51 - 9.02 

3/2 - 5/2 

8912.07 

-0.11 

0.61 

Call 

8912.07 

— 

— 

— 

8927.33 

-0.14 

0.64 

Call 

8927.36 

— 

— 

— 

6587.65 

-0.01 

0.25 

Cl 

6587.61 

5.09eH-06 

8.54 - 10.42 

1 - 1 

7772.00 

-0.23 

0.69 

OI 

7771.94 

3.69eH-07 

9.15 - 10.74 

2-3 

7774.18 

-0.17 

0.63 

OI 

7774.17 

3.69eH-07 

9.15 - 10.74 

2-2 

7775.30 

-0.15 

0.63 

OI 

7775.39 

3.69eH-07 

9.15 - 10.74 

2- 1 

8446.46 

-0.58 

1.65 

OI 

8446.25 

3.22eH-07 

9.52 - 10.99 

1-0 

3905.52 

-0.69 

0.17 

Sil 

3905.52 

1.33eH-07 

1.91-5.08 

0- 1 

4130.88 

-0.12 

0.23 

Sill 

4130.89 

1.74eH-08 

9.84 - 12.84 

5/2 - 7/2 

6347.10 

-0.16 

0.38 

Sill 

6347.10 

5.84eH-07 

8.12- 10.07 

1/2 - 3/2 

6371.37 

-0.10 

0.39 

Sill 

6371.36 

6.80eH-07 

8.12- 10.07 

1/2 - 1/2 

3986.76 

-0.06 

0.20 

Mgl 

3986.75 

7.30eH-06 

4.35 - 7.45 

1 -2 

4057.51 

-0.07 

0.20 

Mgl 

4057.51 

1.02eH-07 

4.35 - 7.40 

1 -2 

4703.01 

-0.06 

0.22 

Mgl 

4702.99 

2.19eH-07 

4.35 - 6.98 

1 -2 

5172.69 

-0.46 

0.33 

Mgl 

5172.68 

3.37eH-07 

2.71 -5.11 

1 - 1 

5183.60 

-0.50 

0.32 

Mgl 

5183.60 

5.61eH-07 

2.72-5.11 

2- 1 

5528.41 

-0.06 

0.24 

Mgl 

5528.40 

1.39eH-07 

4.35 - 6.59 

1 -2 

3895.64 

-0.15 

0.17 

Fel 

3895.66 

9.39eH-06 

0.11 - 3.29 

1-0 

3899.04 

-0.06 

0.26 

Fel 

3899.03 

7.60eH-05 

2.45 - 5.63 

4-4 

3899.70 

-0.22 

0.20 

Fel 

3899.71 

2.58eH-06 

0.09 - 3.27 

2-2 

3900.53 

-0.21 

0.16 

Fel 

3900.51 

7.90eH-06 

3.24 - 6.42 

3 - 3 

3902.94 

-0.23 

0.20 

Fel 

3902.95 

2.14eH-07 

1.56-4.73 

3 - 3 

3903.89 

-0.09 

0.21 

Fel 

3903.90 

7.61eH-06 

2.99-6.17 

4-4 

3906.47 

-0.13 

0.15 

Fel 

3906.48 

8.32eH-05 

0.11 - 3.28 

1 - 1 

3907.93 

-0.07 

0.17 

Fel 

3907.93 

6.67eH-06 

2.76 - 5.93 

3-2 

3917.18 

-0.09 

0.14 

Fel 

3917.18 

4.34eH-05 

0.99-4.15 

2-3 

3920.26 

-0.19 

0.17 

Fel 

3920.26 

2.60eH-06 

0.12-3.28 

0- 1 

3922.91 

-0.32 

0.21 

Fel 

3922.91 

1.08eH-06 

0.05 - 3.21 

3-4 

3927.92 

-0.29 

0.20 

Fel 

3927.92 

2.60eH-06 

0.11 - 3.27 

1 -2 

3930.29 

-0.30 

0.19 

Fel 

3930.30 

1.99eH-06 

0.09 - 3.24 

2-3 

3940.86 

-0.07 

0.21 

Fel 

3940.88 

1.20eH-05 

0.96-4.10 

3-4 

3963.11 

-0.04 

0.25 

Fel 

3963.10 

1.50eH-07 

3.28 - 6.41 

1 -2 

3977.74 

-0.07 

0.19 

Fel 

3977.74 

6.41eH-06 

2.20-5.31 

2-2 

3986.17 

-0.04 

0.18 

Fel 

3986.17 

— 

— 

3-4 

3997.40 

-0.07 

0.17 

Fel 

3997.39 

1.26eH-07 

2.73 - 5.83 

4-5 

3998.04 

-0.06 

0.12 

Fel 

3998.05 

5.70eH-06 

2.69 - 5.79 

5-4 
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Table A.3. Continued. 


^ohs 

(A) 

EW 

(A) 

FWHM 

(A) 

Spec. 

^theo 

(A) 

^ki 

Ej - Ek 
(eV) 

" J/: 

4004.89 

-0.04 

0.23 

Fel 

4004.83 

— 

— 

6-5 

4005.24 

-0.14 

0.19 

Fel 

4005.24 

2.04eH-07 

1.56-4.65 

3-2 

4021.85 

-0.04 

0.12 

Fel 

4021.87 

8.55eH-06 

2.76 - 5.84 

3-4 

4035.65 

-0.04 

0.18 

Fel 

4035.59 

— 

— 

4-5 

4045.81 

-0.28 

0.21 

Fel 

4045.81 

8.62eH-07 

1.48-4.55 

4-4 

4051.93 

-0.05 

0.18 

Fel 

4051.91 

3.30eH-06 

3.40 - 6.46 

3-2 

4062.45 

-0.03 

0.18 

Fel 

4062.44 

1.85eH-07 

2.85 - 5.90 

1 - 1 

4063.60 

-0.20 

0.17 

Fel 

4063.59 

6.65eH-07 

1.56-4.61 

3 - 3 

4071.74 

-0.15 

0.15 

Fel 

4071.74 

7.64eH-07 

1.61-4.65 

2-2 

4107.50 

-0.05 

0.22 

Fel 

4107.49 

1.74eH-07 

2.83 - 5.85 

2- 1 

4132.05 

-0.19 

0.19 

Fel 

4132.06 

1.18eH-07 

1.61 -4.61 

2-3 

4134.67 

-0.06 

0.17 

Fel 

4134.68 

1.25eH-07 

2.83 - 5.83 

2-3 

4143.41 

-0.08 

0.18 

Fel 

4143.41 

2.70eH-07 

3.05 - 6.04 

4-4 

4143.87 

-0.27 

0.23 

Fel 

4143.87 

1.33eH-07 

1.56-4.55 

3-4 

4147.68 

-0.07 

0.29 

Fel 

4147.67 

4.36eH-05 

1.48-4.47 

4-3 

4163.66 

-0.04 

0.15 

Fel 

4163.68 

— 

2.69 - 5.67 

5-5 

4173.46 

-0.25 

0.21 

Fel 

4173.47 

— 

3.01 - 5.99 

2- 1 

4177.66 

-0.12 

0.22 

Fel 

4177.59 

3.71eH-04 

0.91 - 3.88 

4-4 

4191.43 

-0.12 

0.29 

Fel 

4191.43 

2.73eH-07 

2.47 - 5.43 

2- 1 

4202.02 

-0.18 

0.19 

Fel 

4202.03 

8.22eH-06 

1.48-4.43 

4-4 

4216.19 

-0.09 

0.18 

Fel 

4216.18 

1.84eH-04 

0.00 - 2.94 

4-4 

4217.52 

-0.05 

0.25 

Fel 

4217.55 

2.46eH-07 

3.43 - 6.37 

1 -2 

4219.36 

-0.03 

0.13 

Fel 

4219.36 

2.88eH-07 

3.57-6.51 

5-6 

4222.21 

-0.05 

0.22 

Fel 

4222.21 

5.76eH-06 

2.45 - 5.39 

3 - 3 

4227.40 

-0.13 

0.20 

Fel 

4227.43 

5.29eH-07 

3.33 - 6.26 

5-6 

4235.93 

-0.09 

0.21 

Fel 

4235.94 

1.88eH-07 

2.43 - 5.35 

4-4 

4238.80 

-0.05 

0.19 

Fel 

4238.81 

2.41eH-07 

3.40 - 6.32 

3-4 

4247.42 

-0.05 

0.21 

Fel 

4247.43 

1.94eH-07 

3.37 - 6.29 

4-5 

4250.11 

-0.11 

0.21 

Fel 

4250.12 

2.07eH-07 

2.47 - 5.39 

2-3 

4250.79 

-0.18 

0.20 

Fel 

4250.79 

1.02eH-07 

1.56-4.47 

3 - 3 

4260.49 

-0.11 

0.21 

Fel 

4260.47 

3.99eH-07 

2.40-5.31 

5-5 

4271.15 

-0.12 

0.23 

Fel 

4271.15 

1.82eH-07 

2.45 - 5.35 

3-4 

4271.77 

-0.33 

0.22 

Fel 

4271.76 

2.28eH-07 

1.48-4.39 

4-5 

4282.40 

-0.11 

0.20 

Fel 

4282.40 

1.21eH-07 

2.18-5.07 

3-2 

4294.10 

-0.23 

0.20 

Fel 

4294.12 

3.12eH-06 

1.48-4.37 

4-4 

4315.02 

-0.19 

0.28 

Fel 

4315.08 

7.76eH-06 

2.20 - 5.07 

2-2 

4324.99 

-0.04 

0.18 

Fel 

4324.95 

— 

2.18-5.06 

2-3 

4325.77 

-0.16 

0.20 

Fel 

4325.76 

5.16eH-07 

1.61 -4.47 

2-3 

4375.93 

-0.11 

0.18 

Fel 

4375.93 

2.95eH-04 

0.00 - 2.83 

4-5 

4383.54 

-0.27 

0.20 

Fel 

4383.54 

5.00eH-07 

1.48-4.31 

4-5 

4404.75 

-0.25 

0.21 

Fel 

4404.75 

2.75eH-07 

1.56-4.37 

3-4 

4415.11 

-0.16 

0.18 

Fel 

4415.12 

1.19eH-07 

1.61 -4.42 

2-3 

4442.33 

-0.06 

0.21 

Fel 

4442.34 

3.76eH-06 

2.20 - 4.99 

2-2 

4459.12 

-0.05 

0.18 

Fel 

4459.12 

2.52eH-06 

2.18-4.96 

3 - 3 

4461.65 

-0.05 

0.16 

Fel 

4461.65 

2.95eH-04 

0.09 - 2.87 

2-3 

4494.56 

-0.06 

0.18 

Fel 

4494.56 

3.45eH-06 

2.20 - 4.96 

2-3 

4528.61 

-0.11 

0.20 

Fel 

4528.61 

5.44eH-06 

2.18-4.91 

3-4 

4592.05 

-0.05 

0.29 

Fel 

4592.65 

1.61eH-05 

1.56-4.26 

3 - 3 

4920.50 

-0.08 

0.25 

Fel 

4920.50 

3.58eH-07 

2.83 - 5.35 

5-4 

4957.51 

-0.17 

0.45 

Fel 

4957.60 

4.22eH-07 

2.81-5.31 

6-5 

5171.59 

-0.15 

0.30 

Fel 

5171.60 

4.46eH-05 

1.48 - 3.88 

4-4 

5232.95 

-0.07 

0.22 

Fel 

5232.94 

1.94eH-07 

2.94-5.31 

4-5 

5266.54 

-0.03 

0.22 

Fel 

5266.56 

l.lOeH-07 

3.00 - 5.35 

3-4 

5269.53 

-0.27 

0.25 

Fel 

5269.54 

1.27eH-06 

0.86 - 3.21 

5-4 

5270.36 

-0.10 

0.20 

Fel 

5270.36 

3.67eH-06 

1.61 - 3.96 

2- 1 

5328.03 

-0.21 

0.23 

Fel 

5328.04 

1.15eH-06 

0.91 - 3.24 

4-3 

5328.53 

-0.12 

0.29 

Fel 

5328.53 

4.74eH-05 

1.56-3.88 

3 - 3 

5371.48 

-0.15 

0.23 

Fel 

5371.49 

1.05eH-06 

0.96 - 3.27 

3-2 

5397.12 

-0.13 

0.25 

Fel 

5397.13 

2.58eH-05 

0.91 - 3.21 

4-4 

5405.77 

-0.11 

0.21 

Fel 

5405.77 

1.09eH-06 

0.99 - 3.28 

2- 1 
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Table A.3. Continued. 


^ohs 

(A) 

EW 

(A) 

FWHM 

(A) 

Spec. 

^theo 

(A) 

^ki 

Ej - Ek 
(eV) 

" J/: 

5429.70 

-0.14 

0.26 

Fel 

5429.70 

4.27eH-05 

0.96 - 3.24 

3 - 3 

5446.92 

-0.17 

0.29 

Fel 

5446.87 

2.50eH-04 

1.61 - 3.88 

2-3 

5455.59 

-0.09 

0.26 

Fel 

5455.61 

6.05eH-05 

1.01 - 3.28 

1 - 1 

3938.34 

-0.14 

0.23 

Fell 

3938.97 

8.40eH-05 

5.91 - 9.06 

3/2 - 5/2 

4122.61 

-0.16 

0.35 

Fell 

4122.67 

3.30eH-04 

2.58 - 5.59 

5/2 - 5/2 

4178.85 

-0.24 

0.18 

Fell 

4178.86 

1.72eH-05 

2.58 - 5.55 

5/2 - 7/2 

4233.16 

-0.41 

0.19 

Fell 

4233.17 

7.22eH-05 

2.58-5.51 

5/2 - 7/2 

4233.59 

-0.05 

0.19 

Fell 

4233.17 

7.22eH-05 

2.58-5.51 

5/2 - 7/2 

5234.62 

-0.22 

0.25 

Fell 

5234.63 

2.5e-H05 

3.22 - 5.89 

7/2 - 5/2 

4258.16 

-0.10 

0.31 

Fell 

4258.15 

3.10eH-04 

2.70 - 5.62 

3/2 - 3/2 

4273.31 

-0.14 

0.35 

Fell 

4273.33 

9.10eH-04 

2.70 - 5.60 

3/2 - 1/2 

4296.55 

-0.19 

0.19 

Fell 

4296.57 

7.00eH-04 

2.70 - 5.59 

3/2 - 5/2 

4303.17 

-0.25 

0.21 

Fell 

4303.18 

2.20eH-05 

2.70 - 5.58 

3/2 - 3/2 

4351.77 

-0.29 

0.22 

Fell 

4351.77 

4.86eH-05 

2.70 - 5.55 

3/2 - 5/2 

4384.31 

-0.06 

0.18 

Fell 

4384.32 

7.20eH-03 

2.66 - 5.48 

11/2 - 9/2 

4385.38 

-0.22 

0.23 

Fell 

4385.39 

4.50eH-05 

2.78 - 5.60 

1/2 - 1/2 

4416.82 

-0.25 

0.22 

Fell 

4416.83 

2.10eH-05 

2.78 - 5.58 

1/2 - 3/2 

4489.20 

-0.11 

0.21 

Fell 

4489.18 

5.90eH-04 

2.83 - 5.59 

7/2 - 5/2 

4491.39 

-0.15 

0.19 

Fell 

4491.40 

1.89eH-05 

2.86 - 5.62 

3/2 - 3/2 

4508.28 

-0.17 

0.18 

Fell 

4508.29 

7.30eH-05 

2.86 - 5.60 

3/2 - 1/2 

4515.32 

-0.15 

0.18 

Fell 

4515.34 

2.37eH-05 

2.84 - 5.59 

5/2 - 5/2 

4520.21 

-0.12 

0.20 

Fell 

4520.22 

9.80eH-04 

2.81-5.55 

9/2 - 7/2 

4549.51 

-0.31 

0.26 

Fell 

4549.47 

l.OOeH-06 

2.83 - 5.55 

7/2 - 5/2 

4555.88 

-0.21 

0.21 

Fell 

4555.89 

2.26eH-05 

2.83 - 5.55 

7/2 - 7/2 

4576.33 

-0.08 

0.19 

Fell 

4576.34 

6.40eH-04 

2.84 - 5.55 

5/2 - 5/2 

4582.82 

-0.07 

0.26 

Fell 

4582.83 

3.44eH-04 

2.84 - 5.55 

5/2 - 7/2 

4620.52 

-0.06 

0.23 

Fell 

4620.52 

2.53eH-04 

2.83-5.51 

7/2 - 7/2 

4629.33 

-0.13 

0.21 

Fell 

4629.34 

1.72eH-05 

2.81 - 5.48 

9/2 - 9/2 

4666.76 

-0.11 

0.27 

Fell 

4666.76 

1.30eH-04 

2.83 - 5.48 

7/2 - 9/2 

4731.45 

-0.08 

0.23 

Fell 

4731.45 

2.80eH-04 

2.89-5.51 

5/2 - 7/2 

4923.92 

-0.44 

0.29 

Fell 

4923.92 

4.28eH-06 

2.89 - 5.41 

5/2 - 3/2 

5018.44 

-0.51 

0.31 

Fell 

5018.43 

2.00eH-06 

2.89 - 5.36 

5/2 - 5/2 

5169.03 

-0.71 

0.30 

Fell 

5169.03 

4.22eH-06 

2.89 - 5.29 

5/2 - 7/2 

5197.57 

-0.19 

0.21 

Fell 

5197.58 

5.40eH-05 

3.23 - 5.62 

5/2 - 3/2 

5234.62 

-0.22 

0.25 

Fell 

5234.62 

2.50eH-05 

3.22 - 5.59 

7/2 - 5/2 

5264.79 

-0.18 

0.37 

Fell 

5264.81 

3.52eH-04 

3.23 - 5.58 

5/2 - 3/2 

5284.10 

-0.11 

0.23 

Fell 

5284.11 

1.90eH-04 

2.89 - 5.24 

5/2 - 7/2 

5316.65 

-0.34 

0.31 

Fell 

5316.62 

3.89eH-05 

3.15-5.48 

11/2 - 9/2 

5325.54 

-0.07 

0.32 

Fell 

5325.55 

8.00eH-04 

3.22 - 5.55 

7/2 - 7/2 

5425.23 

-0.03 

0.26 

Fell 

5425.26 

9.20eH-03 

3.20 - 5.48 

9/2 - 9/2 

5534.84 

-0.07 

0.20 

Fell 

5534.85 

3.00eH-04 

3.24 - 5.48 

11/2 - 9/2 

5991.35 

-0.02 

0.23 

Fell 

5991.38 

4.20eH-03 

3.15-5.22 

11/2 - 9/2 

6147.69 

-0.03 

0.23 

Fell 

6147.73 

1.30eH-05 

3.89 - 5.90 

3/2 - 1/2 

6238.35 

-0.05 

0.35 

Fell 

6238.38 

7.50eH-04 

3.89-5.88 

3/2 - 3/2 

6247.54 

-0.10 

0.35 

Fell 

6247.56 

1.60eH-05 

3.89-5.88 

5/2 - 3/2 

6456.39 

-0.09 

0.37 

Fell 

6456.38 

1.70eH-05 

3.90 - 5.82 

7/2 - 5/2 

6516.08 

-0.06 

0.39 

Fell 

6516.05 

8.30eH-03 

2.89 - 4.79 

5/2 - 7/2 

3913.47 

-0.23 

0.19 

Till 

3913.46 

1.60eH-07 

1.12-4.28 

7/2 - 7/2 

4028.34 

-0.08 

0.21 

Till 

4028.34 

5.10eH-06 

1.89-4.97 

9/2 - 7/2 

4053.84 

-0.05 

0.18 

Till 

4053.83 

4.20eH-06 

1.89-4.95 

7/2 - 5/2 

4290.23 

-0.23 

0.24 

Till 

4290.23 

4.60eH-06 

1.16-4.05 

3/2 - 5/2 

4300.04 

-0.30 

0.25 

Till 

4300.05 

7.70eH-06 

1.18-4.06 

5/2 - 7/2 

4301.91 

-0.18 

0.25 

Till 

4301.93 

6.20eH-06 

1.16-4.04 

1/2 - 3/2 

4307.88 

-0.27 

0.22 

Till 

4307.90 

4.60eH-06 

1.16-4.04 

3/2 - 3/2 

4312.86 

-0.16 

0.27 

Till 

4312.87 

4.10eH-06 

1.18-4.05 

5/2 - 5/2 

4367.63 

-0.04 

0.19 

Till 

4367.66 

1.90eH-06 

2.59 - 5.43 

7/2 - 9/2 

4395.02 

-0.18 

0.16 

Till 

4395.04 

9.40eH-06 

1.08 - 3.90 

5/2 - 7/2 

4399.76 

-0.09 

0.22 

Till 

4399.77 

3.10eH-06 

1.24-4.05 

3/2 - 5/2 

4417.72 

-0.08 

0.18 

Till 

4417.72 

2.10eH-06 

1.16-3.97 

3/2 - 5/2 
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Table A.3. Continued. 


^ohs 

(A) 

EW 

(A) 

EWHM 

(A) 

Spec. 

^theo 

(A) 

^ki 

Ej - Ei 
(eV) 

" J/: 

4443.80 

-0.10 

0.17 

Till 

4443.80 

l.lOeH-07 

1.08-3.87 

3/2 - 5/2 

4450.49 

-0.10 

0.22 

Till 

4450.49 

2.00eH-06 

1.08-3.87 

5/2 - 5/2 

4468.49 

-0.13 

0.20 

Till 

4468.50 

l.OOeH-07 

1.13 - 3.90 

9/2 - 7/2 

4501.27 

-0.12 

0.20 

Till 

4501.27 

9.80eH-06 

1.12-3.87 

7/2 - 5/2 

4563.77 

-0.07 

0.18 

Till 

4571.98 

1.20eH-07 

1.57-4.28 

9/2 - 7/2 

4571.97 

-0.09 

0.18 

Till 

4563.77 

8.80e-i-06 

1.22-3.94 

1/2 - 3/2 

4589.92 

-0.07 

0.34 

Till 

4589.95 

1.30eH-06 

1.24-3.94 

3/2 - 3/2 

5185.96 

-0.10 

0.29 

Till 

5185.91 

1.40eH-06 

1.89-4.28 

7/2 - 7/2 

3944.01 

-0.30 

0.20 

All 

3944.01 

4.93eH-07 

0.00-3.14 

1/2 - 1/2 

3961.52 

-0.25 

0.20 

All 

3961.52 

9.80eH-07 

0.01 - 3.14 

3/2 - 1/2 

4025.14 

-0.07 

0.18 

CrI 

4025.01 

_ 

_ 

_ 

4077.72 

-0.33 

0.18 

CrI 

4077.68 

— 

— 

— 

4254.33 

-0.22 

0.23 

CrI 

4254.35 

3.15eH-07 

0.00-2.91 

3-4 

4274.80 

-0.20 

0.23 

CrI 

4274.80 

3.07eH-07 

0.00 - 2.90 

3 - 3 

4289.72 

-0.16 

0.21 

CrI 

4289.72 

3.16eH-07 

0.00 - 2.89 

3-2 

4351.77 

-0.29 

0.22 

CrI 

4351.77 

1.20eH-07 

1.03 - 3.88 

4-5 

4541.51 

-0.12 

0.24 

CrI 

4541.51 

3.80eH-06 

3.08-5.81 

3 - 3 

4012.38 

-0.12 

0.17 

CrII 

4012.47 

— 

— 

— 

4558.64 

-0.08 

0.21 

CrII 

4558.66 

8.80e-i-06 

4.07 - 6.79 

9/2 - 7/2 

4588.19 

-0.07 

0.22 

CrII 

4588.20 

1.20eH-07 

4.07 - 6.77 

7/2 - 5/2 

4616.64 

-0.03 

0.24 

CrII 

4616.63 

4.00eH-06 

4.07 - 6.76 

3/2 - 3/2 

4824.11 

-0.08 

0.22 

CrII 

4824.13 

1.70eH-06 

3.87 - 6.44 

9/2 - 9/2 

4848.19 

-0.06 

0.25 

CrII 

4848.23 

2.60eH-06 

3.86 - 6.42 

7/2 - 7/2 

4876.45 

-0.06 

0.26 

CrII 

4876.47 

1.60eH-06 

3.86 - 6.41 

7/2 - 5/2 

5208.47 

-0.10 

0.32 

CrII 

5197.58 

5.06eH-07 

0.94 - 3.32 

2-3 

4030.73 

-0.18 

0.27 

Mnl 

4030.76 

1.70eH-07 

0.00 - 3.08 

5/2 - 7/2 

4033.06 

-0.17 

0.23 

Mnl 

4033.07 

1.65eH-07 

0.00 - 3.07 

5/2 - 5/2 

4034.46 

-0.13 

0.23 

Mnl 

4034.49 

1.58eH-07 

0.00 - 3.07 

5/2 - 3/2 

4110.99 

-0.06 

0.18 

Mnl 

4110.90 

— 

— 

— 


Appendix B: The disk properties of EX Lupi revisited 

Starting with the SED from Sipos et al. (2009), we have fitted the spectral energy distribution (SED) using the RADMC code 
(Dullemond & Dominik 2004). The star is very variable in the optical and IR, and most of the data available are non-simultaneous, 
so the results have to be handled with care. Nevertheless, we can in this way constrain the basic disk parameters, including the 
estimated disk mass. Assuming that the disk is populated by grains with sizes 0.1-100 pm with a collisional distribution, plus some 
optically thin population of small (0.1 to 2 pva) grains in the inner hole (to account for the strong silicate feature), we can obtain 
reasonable fits to the SED with disk masses as low as 1-3x10“^ Mq. The mass is highly dependent on the maximum grain size 
and on the detailed disk structure, which is beyond the scope of this paper. Our estimate of the inner hole in the disk is roughly 
consistent with the values of Sipos et al. 2009 (~0.4-0.5 AU, compared to their estimate of ~0.3 AU), noting that the hole radius is 
also highly dependent on the grain properties and structure of the inner disk rim. In particular, the hole size depends on the sizes of 
the small grains within the hole and whether we assume them to be distributed all over the inner disk or concentrated near the rim. 
Sipos et al. (2009) also demonstrated that a smoothly-curved rim is necessary to reproduce the IR flux. We also find that the best fit 
to the quiescence SED is obtained assuming a stellar radius of 1.5 Rq, slightly lower than the previous value of 1.6 R©, although this 
depends on the quiescence data used, given the small scale variability of the star. Erom this brief analysis, we confirm the presence 
of a dust hole and also suggest that the disk mass does not necessarily needs to be significantly larger than typical disk masses for 
similar CTTS, if we assume a reasonable maximum grain size. 


Appendix C: Accretion rate estimates during quiescence and their effects on the radiai veiocities 

We estimate the accretion rate of EX Lupi using the accretion-related emission lines. Our EEROS/HARPS data are not flux- 
calibrated, and were often taken under non-photometric conditions. Given the variability of the object, it is not even possible to 
estimate the line flux as required to calculate accretion rates (Alcala et al. 2014), so we opt for using archival X-shooter data from 
the ESO database. X-shooter has three spectrograph arms (UVB, VIS, and NIR ), providing simultaneous wavelength coverage 
from ~300 to ~2480nm. EX Lupi was observed as part of program 085.C-0764(A) on May 4, 2010. This epoch is intermediate 
(within a month) between two of our EEROS observations. The slit widths used were Oi'S, Oi'S, and Oi'T, for exposure times of 600 s. 
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Fig. C.l. Line fit parameters for the NC of the Ca IIIR (left) and He II lines (right), colour-coded according to their accretion rate 
(see text). The upper panel shows the radial velocity of the NC line (V^vc), the middle panel shows the amplitude of the NC line 
(A^c). the bottom pannel shows the Gaussian width (cr^c)- There is no evident modulation of the RV with the accretion rate, 
and only the line width for the Ca II NC seems to be higher for the higher accretion rates, although this could be also an effect of 
the presence of a stronger underlying BC. The signatures observed in the NC are thus very robust with respect to variations of the 
accretion rate within one order of magnitude. 


480 s, and 1200 s, respectively, for the UVB, VIS, and NIR arms. We used the ID spectra, which had been calibrated for wavelength 
and flux from ESO Phase 3 Data Release. With the relations between the line luminosity and accretion luminosity, we derive the 
accretion luminosity of EX Lupi using 10 emission lines including He I 5876, 6678, 7065 A, Ha, Call 8498, 8542, 8862 A, Pa/l, 
Pay, and Pacr (Alcala et al. 2014). The mass accretion rates can then be calculated from the inferred accretion luminosities using 
the following relation: 




GM*(1-|^)’ 


(C.l) 


where Rin denotes the truncation radius of the disk, which is taken to be 5 R* (Gullbring et al. 1998). G is the gravitational constant, 
M* is the stellar mass, and R* is the stellar radius. The accretion rates estimated from the 10 lines range from 2 x 10 to 6 x 
10 Mq yr“^ with an average of 4+2 x 10 Mq yr"'. In order to put these values in context with our HARPS/EEROS observations, 
we used the Ha line to infer whether the X-shooter spectrum was taken at a low or a high accretion rate, with respect to the typical 
quiescence values. The Ha 10% velocity width in the X-shooter data is 320 km/s. The two closest EEROS spectra (taken 10 days 
before and 18 days after the X-shooter data) have Ha 10% widths of 320 and 360 km/s, respectively. These values are between 
low and intermediate, comparing the rest of EEROS/HARPS data analysed, although variations in days/weeks are often seen within 
our collection. Comparing the strength of the Ca II BC of the lines, the X-shooter spectrum has also an intermediate-strength BC 
(note the difference in resolving power), also stronger than the two closest EEROS spectra (which have negligible Ca 11 BC). With 
this information, we consider that the measured accretion rate of 4+2 x 10“'° Mq yr“' is probably an intermediate value within the 
typical variations of EX Lupi in quiescence. 

Increases in the Ha EW and velocity wings are usually interpreted as increases in the accretion rate (Natta et al. 2004). In EX 
Lupi, the BC of the Ca II IR triplet becomes stronger when the accretion rate increases. Eor the objects with lower accretion rates, 
according to Ha, the BC of the Ca II IR triplet becomes undetectable. We thus can use the strength of the BC of the Ca II IR lines 
to classify our data according to their accretion levels, considering ’’high accretion” (BC peak > 1 over the normalised continuum 
levels), ’’intermediate accretion” (BC peak 0.2-1 over the normalised continuum levels), and ’’low accretion” (BC peak <0.2 or 
undetectable). These values would roughly correspond to accretion rates of a few times 10“°Mo/yr, between 5xl0“"'-10“°Mo/yr, 
and below 5xl0“*°Mo/yr. This classification allowed us to examine to which extent the NC is affected by variations in the accretion 
rate. Surprisingly, variations in the accretion rate within this range do not seem to affect the NC radial velocities nor amplitudes 
(Eigure [CT| i. Only the width of the NC of the Ca II IR triplet seems to increase in some of the higher accretion spectra, but this could 
be also due to contamination by the stronger BC, especially in the case where the BC is centred around zero velocity. There is also 
no signihcant correlation between the phase and the strength of the accretion rate. In fact, for datasets taken on subsequent days, we 
usually observe that the accretion rate/strength of the features is relatively constant, independently of the shape of the feature. 


Appendix D: Cross-correlation tables for the Ca II line fits 

The following tables include the result of the cross-correlation of the three Gaussian ht for the Ca II 8498 and 8662A lines. Eor 
all the correlations, we perfomed a Spearman rank test providing the correlation coefficient (r) and false-alarm probability (p) for 
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Table D.l. Correlation coefficients for the Gaussian fits to the Ca II 8498 Alines. 


Al 

Cl 

SI 

A2 

C2 

S2 

A3 

C3 

S3 

Al 

0.12/0.45 

0.24/0.12 

-0.08/0.60 

0.31/0.04 

-0.15/0.32 

-0.30/0.05 

0.31/0.04 

-0.35/0.02 

Cl 


0.37/0.01 

0.11/0.48 

0.21/0.16 

0.26/0.08 

0.20/0.20 

0.05/0.76 

0.01/0.96 

SI 



0.41/0.005 

0.50/6e-4 

0.11/0.47 

0.28/0.06 

-0.08/0.61 

0.20/0.18 

A2 




0.56/6e-5 

0.32/0.03 

0.66/2e-5 

-0.38/0.01 

0.59/2e-5 

C2 





0.23/0.12 

0.18/0.24 

-0.13/0.37 

0.14/0.35 

S2 






0.24/0.11 

0.14/0.35 

-0.04/0.78 

A3 







-0.61/9e-6 

0.48/8e-4 

C3 








-0.55/8e-5 


Notes. Autocorrelations in the 3-Gaussian fits for the Call 8498A line. Spearman r correlation coefficient and p probability uncorrelated. Strong 
correlation (false-alarm probability <5e-3) are marked in bold. 


Table D.2. Correlation coefficients for the Gaussian fits to the Call 8662 Alines. 


Al 

Cl 

SI 

A2 

C2 

S2 

A3 

C3 

S3 

Al 

-0.11/0.45 

0.19/0.20 

0.12/0.43 

0.28/0.05 

-0.07/0.66 

-0.14/0.35 

0.07/0.63 

-0.14/0.36 

Cl 


-0.09/0.55 

-0.11/0.46 

0.02/0.88 

-0.02/0.91 

0.01/0.90 

-0.25/0.09 

0.14/0.35 

SI 



0.45/2e-3 

0.21/0.16 

0.40/0.007 

0.48/9e-4 

-0.17/0.25 

0.45/2e-3 

A2 




0.40/6e-3 

0.36/0.02 

0.53/2e-4 

-0.19/0.22 

0.37/0.01 

C2 





0.06/0.70 

-0.2/0.15 

0.09/0.58 

-0.04/0.79 

S2 






0.36/0.01 

-0.11/0.46 

0.1/0.27 

A3 







-0.52/3e-4 

0.68/3e-7 

C3 








-0.63/4e-6 


Notes. Autocorrelations in the 3-Gaussian fits for the 8662A line. Strong correlation (false-alarm probability <5e-3) are marked in bold. 
Table D.3. Cross-correlation coefficients for the Gaussian fits to the Call 8498 and 8662A lines. 


Components 

r/p 

A18498/A18662 

0.86/4e-14 

Cl 8498/C 18662 

0.66/7e-7 

SI 8498/S 18662 

0.63/3e-6 

A28498/A28662 

0.88/2e-15 

C28498/C28662 

0.18/0.23 

S28498/S28662 

0.49/6e-4 

A38498/A38662 

0.81/2e-ll 

C38498/C38662 

0.55/8e-5 

S38498/S38662 

0.54/le-4 


Notes. Cross-correlations in the 3-Gaussian fits for the Call 8498 and 8662A lines. For each pair, the correlation coefficent r and the false-alarm 
probability p are given as r/p. 


the various Gaussian parameters; amplitude (Al, A2, A3), central wavelength (Cl, C2, C3), and width (SI, S2, S3). The numbers 
indicate to which one of the three Gaussians involved in the fit corresponds: narrow component (1), weakest broad component (2), 
strongest broad component (3). For the autocorrelation of every line, we check the relation between all different fit parameters to 
explore potential common dependencies between the three components (Tables [DT| and D.2i. With the two-line correlation, we 


explore whether the independent fit of the two Ca II 8498 and 8662A lines provides consistent results, as both lines of the triplet 
are expected to be similar (Table [d3] i. The most significant autocorrelations appear between the width and amplitude of the lines, 
suggesting that broader lines tend to be also stronger. This is consistent with an accretion scenario where an increased rate requires 
more material coming from more distant regions in the disk and results in the heating of a larger part of the accretion structure. The 
width of the NC could also be affected in case of increased accretion, although this could be also due to contamination by the strong 
BC. In the two-line cross-correlation, we observe that all the parameters are strongly correlated (as expected for the two lines being 
nearly identical) except for the center of the second Gaussian component. This is probably due to the fitting being degenerated; 
the main BC component is often well reproduced by a single Gaussian, so the remaining Gaussian component is often weak and 
poorly defined. Although the three Gaussian model offers very good fits of the line profiles, we observe small deviations from the 
symmetric Gaussian profiles in the BC that suggest a more complex situation. 


Appendix E: The behaviour and properties of individuai iines 

This appendix describes in detail the properties of the most prominent lines observed in EX Lupi. 


The Ca II lines 
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Fig.E.l. Phase-folded curves (using the photospheric absorption line period and phase, which is nearly identical to the He II4686A 
one) for the amplitude of the multi-Gaussian model used to fit the Call IR lines. Top: narrow component. Middle: blueshifted BC. 
Bottom: redshifted BC. There is a marginal amplitude modulation of the BC amplitude. 

The two Ca II IR triplet lines that are covered by FEROS (at 8498 and 8662A, the line at 8542A falls within the 
gap) behave in a very similar way and have similar profiles. They have a very stable NC and a remarkably variable BC, 
which appears redshifted and blueshifted up to velocities over ±200 km/s, and with a peak located up to ±100 km/s. 
The presence of a strong BC depends on the accretion rate, becoming undetectable at low accretion rates. The NC is 
always detected and does not change much (its peak, once extracted from the BC, is always in the range of 1-1.5 over the 
normalised continuum). The lines do not show any significant RV modulation, although there is a clear offset between 
their velocity and the zero velocity of the photosphere. There is also a <1 km/s systematic offset between the velocity of 
the 8498 and 8662A, probably caused by the differences in fhe line opfical depfhs (see main fexf). 

Regarding the amplitude of the BC, there is a low-significance modulafion, especially al limes of high accre- 
fion/strong BC (Figure |E.1| ). The amplilude is maximal al zero velocify, which is the opposite behaviour that we see 
with high-excitation lines such as He II. This is consistent with the Ca II IR lines (and, in particular, the BC) being 
produced in a less embedded area of the accretion column. Therefore, we get more emission when the structure is facing 
us (at zero velocity), while for lines that are produced more deeply in the accretion column, the strongest emission is 
seen when we are looking sidewards to the structure. 

In the IR, we also find fhe Call 8248 and 8927A . Bofh lines are narrow, very weak and fypically noisy, so we cannol 
oblain much information. The 8927A line is not detectable in spectra taken at low accretion. The Call H and K lines in 
the near-UV have a very different morphology compared to the IR lines. The NC±BC structure is not so evident. The BC 
has a very strong blueshifted absorption, which is not smooth, but also appears to have some structure. Some redshifted 
absorption could be also present on certain dates. The blueshifted absorption is clearly correlated with the strength of 
the accretion rate and likely results from an accretion-powered wind. 

Examining the shape of the BC of the Ca II IR triplet vs the phase, we find fhaf they are not completely independent. 
The rapid day-to-day variability observed during outburst (SA12) with strong shifts of the BC peak is also manifest here 
when comparing observations taken on consecutive days (Figure |E!2] ). Selecting only the cases where the BC is strong 
(>0.3 over the continuum level) and has a simple shape, we classified the shape of the peak as ’’centred” (if the velocity 
of the BC peak is in the range ±20 km/s), ’’redshifted” (for velocities >20 km/s), ’’blueshifted” (for velocities <-20 km/s, 
and ’’double-peak” (if the BC has two peaks with height differences <0.2). For each case, we measured bofh the 8498 
and 8662A lines, and the final classificalion depends on fhe averaged resulf, plus a visual inspecfion fo remove uncerfain 
cases (e.g., those where the top of the BC is flat over a large velocity range). Most of the data points are consistent with 
either central or blueshifted BC, given that redshifted components tend to fail the strength criterion, being often weaker. 
Blueshifted BC appear at phases between 0.24-0.5, corresponding to the phases for which the emission lines move from 
no-shift to blueshifted (Figure [£3] ). A KS test shows a probability of 0.9% that the phase correlation is random, which 
points to a significant correlation between the BC shape and the NC line modulation, and thus to a common dynamical 
origin. If we assume that the BC is formed in the same accretion structure where the NC comes from, but at a larger 
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Fig.E.2. Some examples of day-to-day modulation of the BC of the Ca II lines (8498A line in red, 8662A line in black). For each 
set, phase (0) and date (d=JD-2455000) are given. 


distance, the presence of blueshifted BC at these phases would indicate a slightly trailing column, where the material is 
seen as coming to the observer while the footprint of the accretion column where the NC originates moves from zero to 
maximum blueshifted velocity with respect to the observer. 

The He I and He II lines 

We detect several strong He I lines: at 4471, 5015, 5875, 6678, and 7065A. The He I lines differ from the metallic 
lines in that they do not have a clear NC-tBC structure, although this structure was present in the outburst (SA12). They 
are typically broader than the NC of the metallic lines. This was also observed in outburst (SA12) and has been noted 
in other stars (Beristain et al. 1998), with the stronger 5875A line being up to ~20-30 km/s, but not as broad as to be 
regarded as BC. In some cases, a weak BC (which can be blueshifted or redshifted) is also seen. The sharp cut of the 
5875 and 6678A lines on the blue side could be a blueshifted absorption, a possible wind signature. Unlike the Ca II 
lines, the He I lines become sometimes very weak when the accretion rate drops. 

The He I lines show the same kind of modulation as observed in the metallic lines. The amplitude of the modulation 
is of the same order than the RV amplitude for the photospheric absorption lines, although there is a clear offset between 
the zero velocities of both. This is a similar effect to what Gahm et al. (2013) observed for the He I lines of RU Lupi. 
We also observe a mild modulation of the amplitude, similar to what is observed for other metallic lines (Mg I, Fe I, Fe 
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Fig. E.3. Shape of the BC (classified as ’’centred”, ’’blueshifted”, ’’redshifted”, ’’double-peaked”, see text) vs phase. Only data points 
with strong BC and a clear peak are classified. 


II), but with more noise, maybe because the line fits are not so well defined due fo the presence of blueshifted absorption 
components. The potential blueshifted absorption component, which could correspond to a wind, does not show any 
clear RV variation. 

The He II line at 4686A is remarkably strong for a low-mass star. Like the He I lines, it has only one component 
that is broader than the usual metallic NC. The line also appears to have sometimes a redshifted tail (maybe due to an 
infalling, turbulent component), although it could also be caused by a nearby line (there are many weak and strong lines 
in the blue and UV part of the spectrum). He II shows the most remarkable and tight RV modulation pattern of all lines, 
considering that the line is relatively weak and that it is located in a rather noisy part of the spectrum. Its periodogram 
shows a very sharp peak at 7.41d, the same as the period observed for the photospheric absorption lines. There is no 
clear modulation of the line amplitude, although an increase in the blue continuum during episodes of stronger accretion 
may mask accretion-related flux variations. 

The Fe I and Fe II lines 

The Fe II lines are among the strongest lines after the typical H, He, and Ca II ones. In particular the multiplet 42 is 
strong in emission, with good S/N lines at 4923 and 5018 A . The lines show the NC-i-BC structure, although in this case 
the BC is much weaker than for Ca II. The BC shows the same kind of modulation as observed for Ca II. The NC reveals 
a clear RV modulation with an amplitude comparable to that of the photospheric absorption line, but clearly off-phase. 
There is also a small modulafion of the amplitude, of the same type as described for He II (smaller amplitude at zero 
velocity) and thus opposite to the amplitude variations of the Ca IIIR BC. 

The same result is found in other Fe II and Fe I transitions. Other strong iron lines include Fe II at 5362, 5316 A (this 
one is very strong and may have sometimes a weak blueshifted absorption), 4549, 4383, 4385 A and Fe I lines at 5269, 
5270 A. Fe I 5371.49 A has a nearby unidentified line at about 5369.95 A, both of which are weak. The Fe I lines in the 
3700-4000 A region are usually noisy and weak. In general, the behaviour of Fe II and Fe I lines is similar. Their BC 
show bulk velocity shifts that are consistent with those of Ca II. Their NC have clear RV modulations, although some 
of the lines are too weak or too noisy to find a clear paffern in them. Different lines have different zero-point RV and 
amplitudes. 

The Mg I multiplet at 5167, 5172, 5183A 

Several Mg I lines are also detected in the quiescence spectra, among them the multiplet at 5167, 5172, and 5183A. 
The lines are strong and mostly narrow, with BC that are usually faint and noisy but become stronger when the BC of 
Ca II IR gets stronger. They show similar velocity patterns as the Ca II and Fe I, Fe II lines. The NC shows the usual RV 
modulation observed in Fe I and Fe II, opposite to the photospheric lines RV, and the same variations in amplitude vs 
phase, with the lines getting stronger at the maximum blueshifted and redshifted velocity. As for the Ca II lines, there is 
also a small offset in velocity betwen the 5167A line and the other two, although in this case contamination by a nearby 
strong Fe I line cannot be excluded. 
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The OI 7774 triplet and 8446k line 

The O I triplet around 111 A A and 8446 A are also very strong in EX Lupi. The 8446 A complex, being a blended 
triplet itself, has a complex structure, very extended, and rather weak or undetectable at low accretion rates. It does not 
have a clear NC, although it is narrower in the centre and has very extended BC wings (up to ±150 km/s) that changes 
from blueshifted to redshifted following the patterns observed for Ca II. The line could have sometimes a weak redshifted 
absorption, although due to the low S/N and the redshifted emission of the BC, it is not evident. Due to the complexity 
of the line, it is not possible to study its RV. 

The 01 7774 triplet is much stronger and always detected. The 7772 A line is very strong and narrow and easy to fit, 
giving very good results. Having a typical width around 14-16 km/s, it also shows a BC when the accretion rate increases. 
The behaviour of the BC is hard to establish due to the presence of the other triplet lines. Fitting the two weaker lines is 
often a problem due to the contamination by the BC of the strong 7772 A line. The 7772 A line displays the same kind 
of RV modulation as the rest of metallic lines, and has the same weak modulation of the line amplitude. Its width could 
be also modulated, becoming thinner when the amplitude decreases, but this is very uncertain due to the complexity of 
the continuum around the triplet where the BC of all three lines are blended. 


Other lines 

As mentioned in SA12, EX Lupi in quiescence shows more than 200 emission lines, of which at present we have 
only investigated the strongest ones. Other lines include neutral elements (such as Mn I and Cr I), which are weaker 
than the mentioned Fe I and Mg I lines, and other ionized species (mostly Si II and Ti II). Table |A.3' contains the list of 
strong, unblended lines identified in fhe specfra. 

The Cr I 4254 and 4274A lines are relafively sfrong and narrow (~8-10 km/s), allhough Ihey are partially affected 
by pholospheric absorpfion lines. Despile being noisier, Ihey show fhe same RV and amplilude modulalion as fhe olher 
mefallic lines. 

The Si II 6347 and 6371A lines are also relafively sfrong. The 6347A line shows somefimes evidence for a BC wilh 
fhe usual velocily pattern, allhough fhe lines become nearly undeleclable when fhe accrelion rale decreases. Bolh lines 
show fhe RV modulalion observed for fhe resl of mefallic lines, allhough fhe amplilude of Ihe modulalion is slighlly 
lower lhan observed for Ihe pholospheric absorption lines. 

The Ti II lines al 4300, 4302, and 4307A are also detected, despite Ihe poor S/N of Ihe speclra in Ihis area. They are 
narrow (8-12 km/s), allhough one of Ihe speclra shows a hinl of a BC. Despite Ihe high noise level, Ihey show Ihe usual 
RV modulation wilh a peak-lo-peak amplilude aboul 6km/s, and Ihe usual sinusoidal modulation, allhough in blue and 
UV lines Ihe amplilude can be easily distorted by poor continuum determination when Ihe S/N is low. 

Olher lines such as Mn I al 4030A are very fainl, so any allempl to til and exlracl Ihem is noise-dominated. 





